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1. iNTKODUCTiUW 


This report presents the results obtained in an investigation of 
charge-coupled shift register structures. The results which should be 
of use in analyzing CCD designs are computer programs which are 
described in Sections 3, 4, and 5. These programs allow one to obtain 
a two-dimensional electrostatic analysis of the CCD structure to 
determine if the structure allows the efficient transfer of charge. 

Features of significant interest are that two channel oxide thicknesses, 
two levels of metallization, and variable channel doping can be accom- 
modated. The program presented in Section 3 can be used to analyze 
CCD input-output gates , and the program in Sections 4 and 5 can be used 
to analyze 2-phase structures with periodic boundary conditions. 

The remainder of the report deals with the formulation of the analysis 
implemented by the programs. Factors which influence a ccura cy and rate 
of convergence are discussed. 

1.1 Structures of Interest 

Figure 1.1 gives a two-dimensional diagram of a 3-phase structure 
discussed by Boyle and Smith^»^ This structure originally was considered 
as having a uniform channel oxide thickness and uniform channel doping. 
Later, Walden, et. al.^ considered the buried layer in which the doping 
varied from the SiSiO^ interface into the substrate so that the maxi- 
mum potential occurred below the interface. 

Figure 1.2 illustrates a 2-phase structure reported by Krambeck, 
et. al.** in which shallow doping varies longitudinally along the channel. 
Variable doping can preclude problems of charge-trapping by fast-surface 
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Some CCD structures. 


■RfPROBUCIBIIjnY OF THE 
IS POOR 



3 


states by forcing carriers away from the interface and also can preclude 
charge trapping due to potential barriers and wells in the gap and under 
electrode edges. Longitudinal variation of the doping can also give the 
potential asjmvmetry required for unidirectionality of the charge transfer. 
The benefits of these complications of the simple structure are discussed 
in the references cited. 

Figure (1.3) illustrates a 2-phase structure with two-levels of 
metallization and multiple channel oxide thickness. In such a structure, 
overlapping electrodes preclude potential barriers such as may occur in an 
interelectrode gap. The multiple oxide thickness allows the required 
asymmetry for 2-phase operation. An added benefit, is that overlapping 
electrodes shield the channel from static charge which may accumulate on 
the surface of the oxide at the air-oxide interface, thus precluding 
unpredictable variations in the channel potential. Kosonocky and Carnes^ 
have discussed this type of structure. 

1.2 Ideal Mode of Operation 

In this mode of operation the channel is depleted of minority 
carriers except those which are gated into the input end of the channel 
and transferred by the attractive potential produced by the transfer 
electrodes. At the edges of the channel, charge is confined by a 
potential barrier which arises naturally for an N-substrate device 
with a thick field oxide but must be produced by an N+ channel stop 
diffusion for P- substrate devices. The surface potential under an elec- 
trode is obtained from the MOS capacitor formula: 


s 




V. 


FB 





( 1 . 1 ) 
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substrate charge Qg depends upon ipg , the surface potential 
so th^t e<}uution (1«1) becomes a quadratic which Is easily solved for 
i|)g . Q Is the; free-charge, and In the ideal mode of operation, it is 
^mall enough to neglect. A necessary condition for charge transfer is 
that the attractive potential of the succeeding electrodp in the 
direction of transfer must exceed that of the preceding. In the ideal 
operating mode no barriers exist in the interelectrode gap and the 
necessary condition is also sufficient for transfer. 

Equation (1.1) is also useful for estimating the amount of 
charge, Q, which can be accumulated under an electrode before the 
attractive potential is equal to that of the preceding electrode, i.e., 
the full well condition. The amount of charge injected into the channel 
from a junction by means of a control gate can be found in a like 
manner. Extension of the concept by integration allows an estimate of 
trapped charge when potential barriers occur in practical structures. 

In practice, potential wells and barriers may exist and preclude 
complete charge transfer. Their occurrence depends upon the CCD 
structural parameters such as electrode geometry, oxide thickness and 
channel doping. Such phenomena must be investigated by two-dimensional 
analyses. Amelio^ has pointed out the utility of analyzing an assmed 
ideal mode of operation. It is assumed that transfer can occur and 
that Q is sufficiently small to be negligible. Neglecting the minority 
carrier concentration allows the analysis to be carried out without 
simultaneous solutions of the current flow equations. The potential 
distribution calculated is valid for almost a clock period and is a 
steady-state solution provided that charge transfer keeps the channel 
depleted. If the solution shows potential barriers and wells, the 
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results provide a means for estimating the magnitude of the free charge 
which will be trapped in actual operation. 


1.3 Formulation of the Two-Dimensional Electrostatic Analysis 


The complexities which arise in considering a realistic model suggest 
a numerical approach. An approach based on Gauss’ law is chosen because 
it allows the treatment of discontinuities in a straightforward manner. 
There are a number of these discontinuities, including those of the 
dielectric constant, of the conductivity, and those arising from surface 
charge distributions. A rectangular cell structure, such as illustrated 
in Figure 1.4, is most convenient for formulation of a discrete problem 
from a continuous one. The grid points within the cells are chosen to lie 
on the lines of discontinuity. 

A formulation of Gauss* law for each cell should consider a rectangu- 
lar cylinder of unit length, so that: 



D • n dS. 





pdxdy + 




o^dx 


( 1 . 2 ) 

(1.3) 


3-J 

where D is the electric flux density, p is the volume charge density, 
and is a line charge density which always occurs on a line parallel 

with the x-axis. This formulation is consistent with a two-dimensional 
analysis. The electric flux density is defined in terms of the electro- 
static potential u by; 


D = - 


e V u 


(1.4) 
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Figure 1.4 Discrete Cell 


A two point finite difference formula is used to define the gradient 
along each boundary on the cell in all succeeding work. Obviously, 
higher order formulas could be used at the expense of increased compu- 
tation time. Using a two-point formula and keeping in mind that the 
cell point is chosen to lie on any line of discontinuity passing through 
the cell, then one obtains for (1.2) evaluated on the cell: 
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Generally determining requires averaging p over the cell and 

consideration of the fact that p depends in a non-linear manner on the 
potential. Furthermore, o , which will have a specified spatial dis 
trihution, must be averaged over fhe line. These points are considered 
in detail later. 

Equation (1.5) is rewritten in the form given by (1.6) which is 
the typical equation for a nonlinear set. 


"®i,jVl,j '^^i,j'"i,j 


( 1 . 6 ) 


The solution of this set of equations gives an approximation of the 
potential at grid points in the field. It may be noted that for certain 
schemes for finding equation (1.6) is equivalent to the finite 

difference approximation to Poisson’s equation used by Amello^. 


1.4 Survey of Difficulties in Obtaining a Solution 

The major problem is in finding a scheme for solving the large 
system of equations. A secondary, but also Important problem is in finding 
Closely associated with the latter problem is choosing cell sizes 
which give suitable accuracy. 

Let us assume the grid has = I) x points. Then equation (1.6) 
defines a set of equations. Typically this number may range from 

2000 to 10,000 points. Systems of non-linear equations may be solved 
by the Newton-Raphson approach’’ which involves the iterative solution oi 
a set of linear equations in an incremental potential ^^itil all 

6u. .*s are reduced to a suitably small value. However, the resulting 

j 
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set of linear equations is so large typically, that an iterative scheme 
must be used to solve the linear set. Therefore, the scheme may become 
prohibitively expensive if some alternate method will work. This con- 
sideration led us to consider the application of the Gauss-Seidel method® 

t 

directly to the non-linear system. 

The matrix equation may be written as : 


[A] [u] = [Q] (1.7) 

Assume for the moment that [Q] is independent of [u] , which would 
be in fact true except within the semiconductor and on conducting 
boundaries which will be excluded. The Gauss-Seidel method for solving 
such a system is an iterative scheme which gives the (n+1)*^^ approxima- 
tion of Uj as: 


u 


(n+1) 


^ I - 

II J=1 


uCti+l) 
IJ J 


Nx 

— I 

hi j=i+i 


a u^^^ + ^ 
^IJ'^J ^ k. 


( 1 . 8 ) 


where- the single subscript (I) is related to the double subscript 
(l,j) by: 


and. 


I = (j - 1)M + i 


(1.9) 


= M X N 


( 1 . 10 ) 


We do not have assurance that this scheme will converge for a 
non-linear set. (In fact, we have had experiences where it did not 
which are discussed later.)- For a linear set, the scheme can be made con- 
vergent by suitable choice of a relaxation parameter, oj . For w < 1 , 
the method is referred to as under-relaxation and for m > 1 
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over-relaxation. In either case, we have: 

„(n+l) ^ ^(n) ^ ^j;Cn+l) _ ^(n)^ 

where is given by (1.8). The problem here is with convergence 

and obtaining a suitably rapid rate to preclude excessive computing 
time. There is no problem with round-off errors. (See Westlake .) 

Several factors must be considered in the choice of the grid size. 

It can be shown that the differencing scheme used here leads to second 
order accuracy, so that differencing errors are of order 0(h + w ). 

The effect of changing the grid size on the differencing error could he 
easily Investigated experimentally. However, the effect of the grid size 
in determining requires some analytical consideration and is 

treated in detail in Section 2. 
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2, ONE-DIMENSIONAL MODELS OF AN MIS STRUCTURE 


The simple one-dimensional structure is useful for analysis in 
order to determine the effect of grid spacing along the y direction 
in the silicon substrata on the accuracy of numerical solutions. Very 
accurate one-dimensional solutions can be obtained without excessive 
computer time. Furthermore, the physical picture is simpler, and this 
aids in choosing useful models for calculating the charge density. 


2.1 The Depletion-Layer Model 

This model is well known and has proved to be quite useful.^ 

It is illustrated in Figure 2.1, which shows the charge distribution 
assumed. On the silicon side of the Si-Si02 interface, the free 
minority charge is assumed to be accumulated with a surface charge 
density Q. Also present on the oxide-side is the fixed surface sta) 
charge Qgg* Below the surface is a depletion layer extending to a 
depth y = - yjj. The surface normal field in the semiconductor Eg 
and the interface potential ijjg are obtained from Poisson's equation. 


Eg 


1 2 
2 E„ ^d 


The potential in the semiconductor with respect to the bulk is ! 


( 2 . 1 ) 

( 2 . 2 ) 


'J'(y) = j (y - y^)^ 0 5 Y £ + Yd (2.3) 

s 


Applying Gauss's law at the interface to obtain the field in the oxide: 



2.1 Space Charge, Field, and Potential for 
MOS Structure based on Depletion Layer 
Theory. 
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(2.4a) 

(2.4b) 


Integrating to obtain the gate- potential and ignoring the oxiqe cnarge 
and the contact potential of the metal-semiconductor contact, the gate 


po<-fTi*-ial y„ is: 


Vg - (qy d - >5 - 


ox 

ox 


= 


Qb ' Q - Qss 


ox 


+ 


(2.5a) 


(2.5b) 


C A e /e > oxide capacitance, 
ox — ox ox 


(2.5c) 


Equations (2.2) and (2.5) are useful for' making estimates of the 
performance capabilities of MOS structures. 


2.2 Equilibrium Models 

In a more exact solution the hole-electron concentration would be 
allowed to vary continuously according to equilibrium relations s 


p.= 

(2.6a) 

n . 

1 

(2.6b) 


where <j)p and are the quasi-fermi potentials for holes and electrons 

and i|) is the electrostatic potential defined to. be zero at the point 
where p = n = n^ in thermal equilibrium. 

Unfortunately, the current flow and continuity equations must 
generally be simultaneously solved in order to find and (j>^ 
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However, let us consider a transient situation in which a voltage is 
applied to the gate and consider a time interval beginning at a time 
substantially greater than the relaxation time but shorter than the 
excess carrier lifetime. The majority carriers are repelled from the 
interface and minority' carriers are attracted. Since equation (2,2) 
should be a reasonable approximation, we may conclude that the minority 
carriers are exponentially distributed in a very thin region near the 
interface provided the current density' is small. Thus in the substrate 
at a small distance from the interface, Poisson’s equation may be writte 
in terms of nsr 

^ = a. (N e- S' ■ 

dy2 e a 


Normalizing, by kT/q , and y by the extrinsic deBye length 

kTe/q^Ng , we have: 


d^u 

dz^ 



( 2 . 8 ) 


Defining the normalized field by, V = , we may integrate (2.8) to 

obtain ; 


= 2 u + 2(1 ~ e“^) 


(2.9a) 


•’^s dp 

^ a - e-n + 


(2.9b) 


The last term in (2.9a) is the deviation from the depletion layer 

approximation. In most cases of Interest the surface potential will 

satisfy u >> 1 . Suppose u > 5 . The integral in (2.9b) can be 
s s 

broken down into two parts: 
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Z « Z + Z 


(2.10a) 


" /za + 


(2,10b) 


r^s 22 (2.10c) 

" na - e“'^ + 

* 

u 

Figure (2.2) shows the variation of u with respect to Z and also 
the normalized charge density, (1 “ e"’") . Let us check the case 
4- = 1 volt , u = 38.61 . The depletion approximation gives 

S 3 

= 8.787 , which would correspond to AZ “ 3.351 . 
d 

Numerical solutions require consideration of a finite region. 

When we consider applying the preceding models In numerical analysis, 
obviously there will be a significant difference between the ratio of 
the charge to surface potential if we consider interface voltages as low 
as 1 volt. This ratio, r, of the substrate charge for an equilibrium 
distribution to that for a depletion layer approximation is: 


^ (1 + -853 ( 2 . 11 ) 

If one were considering the effects of the surface state charge Qgg 
alone in the absence of an electrode, then significant errors may 
result from the. depletion approximation, particularly if the resultant 
potential is low. However, in this case the absolute error is also 
lowered. This feature seems to save the depletion approximation from 
failure when otherwise it appears grossly in error. 

The space charge density changes from 95% to 5% of Its maximum value 
in approximately 2.3 extrinsic deBye lengths. Therefore, we conclude 
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Chat a space charge^potential' function' model which shows a continuous 
variation will be of significance in a numerical analysis only if the grid 
size, h, is less than an extrinsic de Bye length. This conclusion is 
important because it is typically expensive in terms of grid points to 
maintain the grid size at this small a value along the depletion edge. 

2.3 Choice of Geid Size and Boundary 

The boundary value problem for an MIS structure is most easily 
formulated as one of infinite or semi-infinite extent. Practically, 
substrate thicknesses are much larger than depletion depths. Since a 
finite region must be considered, one must choose a suitable boundary. 

The boundary should be reasonably distant from the region in which the 
potential gradient varies significantly. 

If economy of storage and running time on the computer is para- 
mount, this choice of a boundary should be a subject for careful investi- 
gation. We have chosen the boundary on which the potential Is specified 
(u = 0) at a distance of twice the estimated depletion depth for first 
solutions. By monitoring solutions one confirms that this indeed 
removes the boundary far enough away; however, it shows that a signifi- 
cant amount of the total storage for data points is seldom used. 

The number of points used in the grid can be reduced if one 
chooses a variable grid size. This must be done with care or significant 
errors will result. In the following, we will estimate errors which may 
occur due to round off or quantization of the charge in a cell assuming 
only two charge density levels. Other schemes for averaging the space 
charge may reduce these errors. 
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We choose the boundary at which u = 0 for a maximum surface 
potential, Ug(MAX), at a distance 2Z^(MAX) , where using a depletion 
apprcxitnatioa s 

Z^(MAX) = h Ug(MAX) ■ 

We allow a total of H points to span the distance 2Z^(MAX) spaced 
at (Zj, Z^, •** Using a depletion approximation, we have: 



( 2 . 12 ) 


where Z. is the depletion depth corresponding to any U . We 
d s 

assume that the resolution will be one-half of a cell: 




<2 


n+1 


n- 


,)M 


Therefore : 



\-i-l ~ ^n-1 
2Z 

n 


Let us assume a uniform grid, Z^ - ^n-1 ~ 


(2.13) 


AU 1 

s - _ 


(2.14) 


The percentage error decreases the further the depletion layer extends. 
However, the absolute value of the error increases: 


AU = n h ^ 
s o 

For example, if we choose 50 points to span the grid and estimate that 
the maximum error occurs at n = 25, then we obtain a percentage error of 
about 4%. 
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Let us use (2.13) to define a grid with uniform percentage error of 
100a. This leads to a difference equation: 


with a solution; 


Z - 2aZ„ - Z^ T 
n+1 n n-1 


= 0 


z 


n 


a 


a 


a + *^1 + a^ 
1 + a 


(2.15) 


(2.16a) 

(2.16b) 

(2.16c) 


where in (2.16c) we require an expanding grid and note that a is small 
so that a^ « a. The absolute error is; 


&U 


s 



Thus, it increases exponentially. 

Now let us consider a case of interest where a maximum surface 
potential of approximately 20 volts must be accommodated. Using the 
uniform spacing of 50 points for Z^(MAX) = 40 , we would have 
ho = 2 X 40/50 =1.6 deBye lengths. The maximum error would be 
approximately .8V at the surface with a percentage error of 4%. 

Let us now determine the exponential spacing which preserves the 
same error under the high potential surface. We choose a and N to 
maintain this error to within 1 deBye length of the surface; i. e. , 
ho = 1 . This would require approximately 113 points, of which 94 would 
be within the depletion distance. Changing hp to 2 deBye lengths 
gives 94 points of which 77 would be within the depletion region. Further 
division of the grid within the first exponential step would cost even 
more points. 
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Any non-uniform spacing which preserves the accuracy for high 
potential at lower potentials will require more points in the grid. The 
grid points may be better utilized (i.e., fewer zeros stored) with the 
non-uniform grid. However, there is another consideration which weighs 
heavily against compromising the accuracy of solutions for higher sur- 
face potential to obtain accuracy for lower potentials. 

In analyzing a CCD structure, the differences in potential along 
the interface will be of great interest. It will be essential to pre- 
serve accuracy in finding the potential in the high potential region; 
otherwise, a relatively small percentage error in this region may lead 
to a relatively larger error in the potential gradient along the surface 
than the same percentage error at a region of lower potential. 

The conflicting requirements for accuracy when the depletion region 
is both deep and shallow and the need to economize the number of grid 
points led us to a compromise. The larger region of the grid is uniformly 
spaced, extending a distance of K • > where the first run value for K 

is 2. Subsequent runs may reduce K if too large a portion of the field 
is unused. This uniform grid will have cell height typically between 
1.5 and 3 deBye lengths, although by some experimentation it can sometimes 
be reduced to less than a deBye length without an excessive number of 
points. This uniform grid is spliced into an exponential grid extending 
from the surface to a distance of deBye lengths, defined by: 


h 

n 


h 

0 


n-1 

cc 


Z 


n 


n 

k=l 


n 


- h 

o^ a - 1 


(2.17a) 


(2.17b) 
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Choosing, a -f 1 + a , where a is small, we obtain a fractional error 
which is 1 for n = 1 and approaches the value "a" for large n . For 
the examples which we have considered with doping between 10^^ and 
5 X 10^^, oxides from .1 to .5 ji, and voltages from -1 to 20 volts, we 
have observed that the depletion layer is usually at least 10 deBye 
lengths from the surface. The uniform grid usually insures maximtun 
errors in the range of 10 w 30%. .Now in order to make a smooth tran- 
sltion across the Interface, i. e. , no large changes in coefficients of 
the equations for , we require hp to take certain values related 
to the grid size in the oxide. These values, in deBye lengths, range 
from .375 (low doping) to 2 (high doping). In the former case approxi- 
mately 37 points are required to extend 10 deBye lengths and obtain 10% 
error bound. In the latter case, these requirements are Incompatible 
and the exponential grid is not used. Thus far we have used an 
exponential grid which extends 2 deBye lengths. 

We have discussed in detail the effect of grid size on round-off 
error, assuming that binary values can be assigned in each cell half. 

The grid size of course affects the solution otherwise, depending upon 
the order of the differencing schemes used. We have used the previous 
considerations to aid in choosing the grid scheme, and then we have 
experimented with the grid size, halving and quartering it, to determine 
the effect upon convergent solutions. When these solutions are in 
agreement, we feel the techniques produce numerical results at least as 
good as the model, and, consequently, of use in determining design 
feasibility. 
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3. TWO-DIMENSIONAL ANALYSIS OF 3-ELECTRODE STRUCTURE 


The 3-electrode structure is illustrated in Figure 3.1, and it could 
be considered as a storage cell for a 3-phase CCD structure or as a 
gating structure for a CCD or SCT (surface charge transistor). A zero 
potential boundary encloses the structure to obtain a finite region 
for analysis. The effect of the location of this boundary can be 
determined empirically. The program given in Appendix A for analyzing 
this structure is based on the Gauss-Seidel iteration procedure modified 
by use of a relaxation parameter. The programmed formulas are: 
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(3.3c) 

C. 
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(3.3d) 


The iterative procedure starts one line from the bottom boundary and one 
line to the right of the left boundary. It continues to the right along 
a row, then back to the left -and upward. All points are considered 
except those at which' the potential is fixed. After the field is swept 
through, the procedure is repeated until satisfactory convergence is 
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3.1 Three-Electrode Structure, 
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obtained. The coefficients in (3.3) are in fact calculated from the 
relative dielectric constants and is divided by the free 

space permittivity. 

We have not used a variable x-grid spacing although this is simple 
to accomplish. This can be done easily by modifying the subroutines 
COEFF and RELAX. However, it will require the operation of COEFF at 
each point rather than once for each line. This will considerably 
increase computation time, and, in our opinion, unnecessarily so for 
most purposes. 

The cell charge in the substrate, , is calculated using the 

subroutine AVERO. Several algorithms were evaluated and will be dis- 
cussed in subsection 3.2 


3.1 Grid Spacing 

As mentioned above a uniform x-grid was used with various values 
for w evaluated. Figure 3.2 Illustrates the y-grid spacing parameters. 
Spacing, through the electrodes and oxide are uniform with ratios 
roughly equal to the ratio of the dielectric constants,. Above the 
electrodes the spacings are exponentially increased, with the spacing 
through the electrodes serving as the exponential base. 

The spacing through the oxide serves as the base for the spacing 
in the substrate immediately below the interface. The basic step size 
is : 


h 

o 


(t /L3)(e^/e ) . 

ox s ox 


Steps are exponentially spaced for a distance of 2 or more extrinsic 
deBye lengths. These exponential steps are calculated by the subroutine 
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YSPAC. The objectives in choosing this spacing were two-fold. First, 
it is desirable to match the coefficients and , which are large, 
on the boundary of the Si-Si02 interface. Second, we would like to 
Increase the accuracy of computing on the depletion edge for 

shallow depletion when it is possible to do so without compromising 
the accuracy on a deep depletion edge. The spacing below this region 
is uniform because uniform spacing spans the ’largest distance with 
the smallest maximum spacing, an obvious but important point when 
considering calculation of on a deep depletion edge. The program 

sets phe cell size equal to a constant, whose default value is 2, times 
the maximum depletion depth (estimated) divided by an integer (input 
data). We have used values of 25 and 50 for the integer. 

3.2 Calculation of Q. . 

Several algorithms were tried for AVERO. The first one utilized 
the exponential dependence on u of p , the majority carrier approximation, 
discussed in Section (2). The average potential in each quarter of a 
cell was used, and the sum of the charges in the four quarters obtained. 
The average of top and bottom halves were used similarly, and then, 
finally, the average of the cell used. 

There were difficulties with all these methods- First, when many 
iterations are required for convergence, l,e,, starting with a poof 
approximation, oscillations of the potential in the substrate occur. 

Once the potential changes sign, large change in the exponential may 
occur. Reduction of the relaxation parameter damps this oscillation 
but slows convergence. If the exponential is made to depend on the 
negative of the potential magnitude, a new type of phenomena occurs which 
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is worse. The depletion edge continues to advance with a negative dip 
in the potential (for a p-substrate) occurring. A stable solution will 
finally be obtained, which is obviously useless, where all the substrate 
is depleted. 

The cure for the problem is to preclude the substrate potential 
from taking the wrong sign. When this step is taken, we find that the 
results obtained with an exponential function does not, in most cases, 
warrant the increased computing time. However, as pointed out earlier 
in section (2) , it must make a difference when the depletion depth is 
shallow compared with an extrinsic deBye length. In most cases, we have 
used a simpler averaging technique .involving binary values for the top 
and bottom of the cell. Whether binary or exponential weights are used, 
it makes little difference whether averaging is done on the right and left 
sides of the cell. 

When the interface is considered, then the surface state charge 

Q must be included as well as any charge due to shallow doping. The 
s s 

shallow doping varies along the channel, and averaging on the right and 
left hand sides of the cell should be done. 

3 . 3 Convergence 

In the following we shall refer to the residual, R, which is the 

change in the potential at a point during an iteration. is the 

maximum change on the interface and R is the maximum change on the 

max ^ 

entire field. As a test case we chose a configuration similar to 
Amelios® with the following parameters: 




3.2 Interface Residuals vs. 
Iteration Number. 
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3.4 Surface Potential Profile for Various Number of 
Iteration. Label (b) indicates starting 
solution with oxide field equal to zero. 
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We concluded that it would take about 240 iterations, using 193 
seconds of CPU time on a UNIVAC' 1106, to develop the essential features 
of a solution using 8736 grid points and a zero initial estimate for the 
potential. Starting with a better initial estimate the time is cut in 
half. Using a coarser grid speeds the convergence with the penalty of 
loss of resolution. Figure 3.6 shows the solution in the vicinity of 
the two electrodes with higher voltages. Reduction of the y grid from 
91 lines to 53 lines shows the gap potential slightly higher but 
essentially the same solution with about one-half the computing time. 

A reduction of the x-grld from 96 to 48 lines shows essentially the same 
solution. This latter solution uses 2544 grid points and the CPU time 
is 80 seconds. The maximum y grid spacing in the substrate is 0.49v com- 
pared with a deBye length of 0.183!j. However, the agreement of all the 
solutions under the low voltage electrode indicates that the space charge 
averaging technique is adequate. 

3.4 Program Features 

The program allows the operator to establish a large number of the 

significant parameters of a CCD structure. Geometrical features such as 

oxide thickness, electrode thickness and width, and gap spacing may be 

controlled. Other parameters such as Q , substrate doping, shallow 

s s> 

channel doping (treated as a surface distribution) , and the electrode 
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voltages can be established. The number of steps through the coarser 
substrate grid, the oxide, through the electrodes and in the air space 
above is under the operator's control. The number of steps through 
the fine substrate grid is set in the program. These parameters are 
read in as input data. 

The output of the program gives the estimated maximum depletion 
depth, the deBye length, the y-grid values with respect to a zero origin 
at the interface, the indices of boundaries of discontinuity, the 
residuals and , and potentials along chosen lines. 

The number of iterations is fixed by the parameter ITER, and the 
iteration numbers at which a printout of u is desired is fixed by the 
parameter LAP. The nxmiber of lines to be skipped is fixed by the 
parameter JC. These parameters are read in as data. Anyone with pro- 
gramming experience can easily modify these features to suit his own 
ideas of convenience. We believe that output of such a program should 
be monitored to determine the nature of convergence rather than to rely 
on built-in tests. Simple built-in tests would be easy to implement 
but may not be reliable or more complex tests may result in unnecessary 
costs of computing time. 

Instructions on preparation of data cards for using the program 
are given in Appendix A. 

3.5 Application of the Program 

The applications of the program to 3-phase structures are obvious. 
For example, the potential profiles can be checked for wells and 
barriers as a function of gap width and channel doping. Let us consider 
the structure for which results were given in Figures 3.5 and 3.6 and 









33 


assume a practical clockwave form for the transfer electrodes. Such a 
waveform is shown in Figure 3.7 which shows a rapidly rising leading 
edge and a slow trailing edge. By analyzing the potential profiles for 
several voltage values on electrode 2, we can find the range in which 
transfer can take place and estimate the fraction of a clock period for 
which transfer can occur. Profiles are shown in Figure 3.8. These 
solutions are obtained using a 2448 point grid which preserves the 
essential features of the solution as was illustrated in Figure 3.6. 

The solution for = 4 used 240 iterations, although inspection of 
the output data showed that 120 iterations gave the same answer to within 
0.6%. Succeeding values for used 80 iterations each. The dotted 

curve shows the solution starting with V 2 = 10 and running 240 itera- 
tions. The entire set of data can be generated with 182 seconds of CPU 
time. The results show that charge can be transferred during the time 
is between about 5 and 7.5 volts. 

The program is also applicable to the study of a 2— phase structure, 
figure 3.9 shows the potential profile at the interface for a structure 
with doping in the gap and under the left hand electrode edges such as 
proposed by Krambeck et. al. (Section 1, Ref. 4), The results show how 
that doping effects an asymmetry which allows unidirectional transfer of 
charge. The doping and geometric parameters were established on the 
basis of one-dimensional calculations, and the program verifies that the 
scheme will work. In this case, the experimental work at Bell Labs has 
proven the feasibility. 

The last application considered here is for studying gate action. 

We give no results, since it requires modification of the program. We 
plan to add the modification later to be used as an option, but a 
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3.7 Profiles for Implanted 2-Phase Structure. 
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programmer can do it for himself. Suppose that a diffused junction is 
to the left and under electrode 1. If the junction is forward biased, 
the potential in that region is very near the substrate (about - .5 volts). 
The potential can be set to zero in that region, and the effect of varying , 
the voltage can be studied. On the right hand side a collector 

junction region can be located where a positive (for n— channel) potential 
is specified, and the effect of V 3 in controlling the output gate can 
be studied. The techniques used in specifying the gate voltage and in 
carrying out the relaxation through the electrode region are 
applicable. 

3,6 Determination of the Electric Field 

We have not Included an algorithm in the program for calculating 
the electric field. In the case where one wishes to estimate transit 
time, one must, of course, have this data. It is on this point that some 
criticize the numerical solution method because simple two-point 
numerical differentiation leads to large discontinuities in the field 
profile. The method for estimating the transit time must be chosen 
before it is clear just what data is needed; however, for the preceding 
cases only the x-component of the field at the interface will probably 
be sufficient. After a solution for Ug is determined, a high-order 
polynomial fit to the potential profile will allow a smooth approxi- 
mation of the field at grid points. 
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4. TWO-DIMENSIONAL ANALYSIS OF A TWO-ELECTRODE STRUCTURE 

The two-electrode structure is a natural one for studying charge 
transfer in a 2-phase CCD. The program given in Appendix B applies the 
methods discussed in Section 2 to analyze the structure shown in Figure 
4.1. The program has a few features which are different insofar as 
the programming is concerned and allows the consideration of an oxide 
with two thicknesses. Basically the progr^ is the same as tnat in 
Appendix A. Periodic boundary conditions are natural for studying the 
inter and intra-cell transfer of charge, and these are incorporated in 
the program. 

4 . 1 Program Features 

The two oxide thickness feature mentioned above is a completely new 
and distinct feature of this program compared with the 3-electrode pro- 
gram. The width of the thin oxide region is controlled by the parameter 
LOXl. The location of the thin oxide region with respect to the left 
hand edge of an electrode is controlled by the parameter LSPOX. The 
thin and thick oxide thicknesses are TOXl and T0X2 respectively. 

In addition, the shallow doping feature is slightly different from 
the 3-electrode program. In the 3-electrode program the channel charge 
was lumped with QSS to form an effective interface surface charge. In 
the 2-electrode program QSS and QS(I) are separated. The doping is 
Included in QS(I) and is read in as an effective surface charge, QP under 
the pads and QG in the gaps. These numbers are positive if the ionic 
charge is positive. Internal to the program, and beginning with line 194 
of the program listing after the instruction labeled 126, is an algorithm 



4.1 Two-Electrode Structure, 
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which determines a voltmie distribution coefficient CHION(J). The algo- 
rithm distributes the ions over 10 lines of the y-grid according to a 
Gaussian formula with a mean of YEAR and standard deviation of SIGMA. We 
realize that this is a simplistic procedure, especially if two types of 
ion are involved. However, this procedure seems to obtain solutions 
with essentially correct features without introducing kinks in the 
potential profile due to abrupt changes in QS(I). The product QS (I)*CHI0N(J) 
gives the volume space charge weight in a cell. This technique can be used 
to handle buried layers where the doping is of a polarity and sufficiently 
heavy to shift the maximum potential away from the interface. A trial 
run may be required to determine whether more than 10 lines are needed 
for the distribution. If the layer is deep compared with a deBye length, 
which doesn't seem likely, then shifting the desired 10 (or whatever) 
lines to a location beginning below the interface may be desirable. This 
can be done several ways in the program. 

The doping QP starts LSPAC steps from the left hand electrode edge 
(LSPAC negative places QP to the left) and extends for a width of LDOPE 
steps. Making LDOPE = LP and LSPAC = 0 and QP = QG, or simply letting 
QP = QG = 0 and QU equal the uniform implanted layer surface density, 
allows treatment of a buried layer. 

The 2 electrode formulation allows a finer grid for the same electrode 
and gap geometry, and the program has a feature which allows the relaxa- 
tion procedure to bypass the lower lines in the substrate where changes 
are extremely slow. As changes occur which advance the depletion layer, 
the relaxation procedure then start a line lower, etc. All of these 
factors allow the program to run with a higher resolution grid and shorter 
time when compared with the 3~electrode program. 
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4.2 Application of the Program 

We have chosen an example which we believe illustrates the manner 
in which two-dimensional analysis programs are useful as design aids. 
Creative thinking of a semi-quantitative nature concerning CCD's will 
probably rely on one- dimensional models. The two-dimensional analysis 
can check the validity of the conclusions arrived at from a one- 
dimens ional cons Ider ation . 

The structure shown in Figure 4.2 appears to have some promise for 
2-phase CCD application. The two-oxide thicknesses allow an asymmetry 
necessary for unidirectional transfer. The electrode voltage will pull 
the surface potential higher under the thin oxide. The surface charge 
in the gap region, from either QSS or uniform implantation, pulls the 
surface potential high enough to preclude a barrier in the gap. Implan- 
tation through the thin oxide allows a sufficient surface charge to 
maintain a high potential on the right side of the interface even when 
the electrode voltage is low. 

We start with the MIS equations : 
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where in (4.2) it is assumed that the oxide field is zero. Now we 
assume that Qs^ ” Qg > ^sl ~ "^^g ^ ^oxl ” ^ 0 x 2 * Next we define 
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the parameters g = qNeg/CQ^ 2 Qg then divide ijjg 

by i^»g and write expressions for the normalized potential in each 
region : 


* . ^ . o+ll j + 2gi.26(<.^.2S^£ 1/2 
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Let us choose the oxide thicknesses to be O.ly and 0.5 m , thin enough and 
thick enough for selective ion implantation. Now with zero electrode 
voltage we wish to make so that minority carriers are 

shifted always to the right hand edge. Setting a = 0 and solving for B 
from (4.4) and Y from (4.3), we obtain: 
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We choose the positive signs in equations (4.5) and (4.6) and let 

“ 3-/2 , = 3/4 . This gives B = 0.853 and Y = 3.066 so that 

Q /q = 6.7 X lO^^cm ^ ^ = 8.6V , and Q /q = 2.05 x lO^^cm ^ with 

6 8 P 

a substrate doping of N = 3.8 x lO^^cm ^ and n = 5 . 

Figure 4.2 shows the normalized surface potential vs. the normalized 
transfer electrode voltage. Choosing V = 20V gives a = 2.33 , and 
this should work well for obtaining a transfer. In Figure 4.3 the 
maximum potential profile is shown along with the location of the maximum 
potential. In the more heavily doped regions the maximum potential occurs 
below the interface as is shown in Figure 4.3. This accounts in part for 



4.2 Normalized Potentials on Right and Left Hand Sides of 
Electrode Versus Normalized Gate Voltage. 
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the rise of the potential above the predicted value from the one- 
dimensional approximation shown in dotted lines. A gaussian distribu- 
tion of the implanted charge is assumed. It is observed that a potential 
barrier occurs at the right hand edge of the low voltage electrode for 
V = 0- Therefore, charge would be trapped under this electrode. One 
observes that the potential maximum is located deeper underneath the 
low voltage electrode than under the high; therefore, the ratio of the 
trapped charged to the amount which can be stored under the high voltage 
electrode is somewhat lower than it first appears. Most of the charge 
which can be stored under the high voltage electrode will be on the 
right-hand side with the higher oxide capacitance. Neglecting the 
differences in location of the potential maximum the ratio of the trapped 
charge to maximum stored with electrode voltages of 0 and 20V is : 


2.3 C , 

0 x 1 

5.2 C , + 15.2 C , 

0X2 0X1 


= 14%- 


Figure 4.4 shows the effect of misalignment of the electrodes with 
respect to the oxide: Misalignment to the right accentuates the barrier, 

and misalignment to the left introduces a well for trapping charge. The 
absence of a peak in the profile at the right hand edge of the 20V elec- 
trode was unexpected and caused us to do some intensive checking of the 
program but we found no errors. Finally, Figure 4.5 shows the profile 
for a sequence of values for indicating the possibilities for 

charge transfer for practical clock waveforms on the transfer electrodes. 
In conclusion, even with + In misalignment the trapped charge is of a 
tolerable magnitude for CCD operation. . 
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4.3 Possible Modifications 

Some may have objections to the using the smaller increments 
beneath the interface. This feature can be modified by rewriting the 
subroutine YSPAC. The feature can be deactivated by removing the card 
sequence; HO = 3.0*T0Xl/L3 to J2 = Jl + L2 and setting J2 == Jl. 

The distribution of implanted charge can be modified by allowing it 
to go further into the grid, i.e., increasing the limit LL = 10 and 
dimensioning CHION( ) suitably. The function used can be changed by 
defining parameters other than YEAR and SIGMA or changing their 
values. The most flexible approach is to substitute a subroutine function 
say FIMPL(A,B, for the gaussian function used. However, a 

suitable function can be incorporated in the main program also. 

4.4 Use of the Program 

The program uses approximately 20,000 words or yu,uuu BYTES of 
storage. Of this 11,000 words are used for storage of the field points. 
This can be reduced by a factor of 4 for many cases by redimensioning 
the arrays for a coarser grid. The run time for 240 iterations and a 
3504 point field is 240 seconds of CPU time. 

The charge entered as QG, QP, AND QS will be distributed 
according to the same function as the program is now written. QG will 
be located in the gaps between the electrodes. QP can be placed any- 
where under or to the side of an electrode. QU is uniform all across 
the surface. QSS is treated as a pure surface charge uniformly 
distributed. 

The thin oxide region can be located as desired. The left hand 
edge is specified and the width of the region specified. 
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Read in data for the program is as follows : 


I; READ(S,1) JA,JB,JC,MAXU,IRES 

1 FORMAT (SIIO) 

JA: lowest line of potentials to be printed out. 

JB: highest line of potentials to be printed out. 

JC: steps between lines printed out. 

Internally, program adjusts JA^l , JB_<J6 , such that the line 
for J=J2 (interface) is printed. 

MAXU : set to integer greater than zero if the maximum 

substrate potential profile is desired. 

IRES; set to integer greater than zero if a printout of the 
residuals after each iteration is desired. 

(Recommended for a first time through until 
confidence in convergence is established.) 

2. READ(S,2)T0X1, T0X2,TEL,W,HFAC 

2 FORMAT (5F10. 3) 

TOXl: thin oxide thickness in microns. 

T0X2: thick oxide thickness in microns. 

TEL: electrode thickness in microns. 

W: width of grid cell in microns. 

HFAC; (HFAC > 1.0) multiplies estimated depletion 
depth to set zero potential boundary in substrate. 

Default value is HFAC = 2.0. 

3. READ(5,4)QSS,CSUB,QG,QP,QU 
4 FORMAT (5E10. 3) 

QSS: surface state charge, cm (i.e. QSS/q) 

O 

CSUB: substrate doping, cm” . 

QG; gap doping, cm“2. 

QP : under electrode doping , cm 

QU; uniform doping, cm”^. 

4. READ(5,6)V1,V2,ITER,LAP 

6 FORMAT (2F10. 3, 2110) 

Vl: left hand electrode voltage, volts. 

V2 : right hand electrode voltage, volts. 

ITER: total number of iterations. (Try 120-240 for start)- 

LAP: printout control, printout occurs when Iteration 

parameter LOOP /LAP = Integer 

5 . READ (5 , 1 8 ) LP , LG , LOOPE , L SP AC , LOXl , LSPOX 

LP: Pad width in units of w. 

LG: gap width in units of w. 

LDODE; width of underpad doping in units of w. 

LSPAC: location of underpad doping from left hand 

electrode edge in units of w, positive to the right. 

LOXl: width of thin oxide region in units of w. 

LSPOX; spacing of thin oxide region with respect to left 
hand electrode edge in units of w, positive to the right. 
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6. READ(5,408)OM1,OM2,OM3,OM4,OM5 
408 FORMAT (5F10. 3) 

These are the relaxation parameters for the 5 regions. After 
early experimentation we used: (1.8, 1,8, 1.5, 1,5, 1.5). 

The user may start with these and do his own experimentation. 

7. READ(5,6) V1,V2,ITER,LAP 

Repeat for as many subsequent values as required. 

The solution builds on the last solution obtained. 

If VI and V2 are not changed more than IV 
ITER = 50 should suffice. 

8. BLANK CARD. 

Computation terminates after test on a blank card. 
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5. TWO-DIMENSIONAL ANALYSIS OF A FOUR-ELECTRODE STRUCTURE 

This structure, illustrated in Figure 5.1, seems to be popular 
for CCD design because it provides complete electrostatic shielding of 
the channel. Originally, it was recognized as a possible solution 
to the problem of electrostatic barriers which may occur in the gaps ; 
however, it was later recognized that shielding- precludes the build-up 
of a static surface charge which influences the interface potential in an 
unpredictable manner. 

Although we have seen no reports of operation of CCD's constructed 
in silicon-on-sapphire (SOS) , we believe that this type of design is 
under consideration. Imaging devices using a 4-electrode structure would 
appear to be quite attractive as an application of silicon on transparent 
sapphire. SOS has not been attractive for bipolar devices because of the 
relatively low life-times due to surface recombination. Apparently lack 
of availability of the material or economic factors have held back 
application in the MOS area.. We did a rough sensitivity check using the 
surface state density data of Elliot and Anderson^ and assuming a 2.5p sec 
lifetime and that the thin silicon layer would be depleted. For a one 
micron layer, visual Inspection reveals that most of the incident light 
is absorbed; therefore, we conclude that after a reflection from a con- 
ducting boundary essentially all incident radiation would be absorbed. 

We obtained a sensitivity of approximately 4mw/cm^ assuming a quantum 
efficiency of unity. After making several conversions of units which we 
hope were correct, we found that the sensitivity of a "sensitive" phot- 
transistor currently available is approximately 2mw/cm^. This latter 
figure no doubt includes the effects of a lens, etc., but is still useful. 
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5.1 Four-Electrode Structure, Silicon is assumed to 
be depleted by an appropriate arrangement not 
revealed by two-dimensional geometry used for 
periodic boundary conditions. 
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We assume that others have checked the sensitivity and found it to be 
tolerable. 

The program which we describe in this section can be used to 
analyze a 4-electrode SOS structure. However, it can also be used to 
analyze a 4-electrode CCD structure on a silicon substrate. 

5.1 Program Features 

The program provides for a sapphire substrate. Since the substrate 
is relatively thick compared with the epitaxial silicon layer, the y-grid 
spacing in the substrate is exponential. Treating the sapphire as 
infinite in extent and applying LaPlace's equation we estimate that the 
field is attenuated by e at a depth equal to the width of one cell 
in a periodic structure. Periodic boundary conditions are used for the 
potential variation along the x-directlon. Exponential y-grid spacing 
would correspond roughly to equi-potential points along an x-cut through 
the grid. The substrate charge gives rise to a "d.c." field component, 
the value of which depends upon the thickness of the sapphire substrate 
and the boundary condition of the surface. It seems probable that the 
sapphire would be greater than 250p thick and that a transparent, 
grounded conductor would be placed on the- surface to eliminate static 
charge. With this thickness the d.c. field component would be small; 
therefore, we set the field equal to zero at a distance of one cell 
depth into the sapphire. For the silicon-substrate option, the boundary 
conditions in the substrate are similar to the previous programs. 

The analysis assumes that the silicon is depleted; therefore, the 
results are meaningful only if depletion can be maintained. It is of 
course conceivable to have a silicon film thick enough so that depletion 
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is not complete. Such a structure may be in fact required to allow 
depletion for the full length of the channel, but we do not know the 
answer at this time. If this is true, then the sapphire substrate 
feature is superfluous; however, the program can still handle the prob- 
lem. Since total absorption is possible within a thin, completely 
depleted film, it seems desirable to obtain this situation if possible. 
The program allows the dimensions of the electrodes, the doping, oxide 
thickness, electrode thickness and electrode voltages to be controlled. 
The relative positions of the electrodes is fixed; however, this seems 
to be. no serious limitation. The number of iterations is controlled on 
the same data card with the electrode voltages so that subsequent runs 
with perturbations of the electrode voltages are possible. 

Spacing is uniform in each of six regions and exponential in the 
sapphire and airspace regions. The exponent will usually be smaller for 
the silicon substrate, since the grid usually will span a smaller 
distance. In this case, the maximum error which may occur due to round- 
off of the space-charge is; 



.2 - 1 ) 
L2(a-1) + - 1 


If one is concerned with error other than at the interface, the above 
relationship is not valid. We feel that the grid spacing for the 
silicon substrate needs further consideration, since it is definitely 
compromised for simplicity and compatibility with the case for a 
sapphire substrate. 
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5.2 Application for the Prograin 

Figure 5.2 shows results obtained in applying the prograin to an SOS 
structure. Curves are shown for zero space charge, for zero pad volts, 
and for the composite. Since the problem is linear, the two individual 
source curves should add to produce the composite. It can be seen that 
they do. One notes that the peaks of the potential curves for zero volts 
are not equal as they should be. It was suspected and confirmed that a 
programming error resulted in shortening the upper left electrode and 
shifting of the upper electrodes to the left. This gave complete 
coverage for the center electrode but not the left hand and right hand 
gaps. The error was corrected. The results, however, do not give the 
asymmetry required for transferring charge undirectionally. 

Figure 5.3 shows the potential profile for electrode voltages which 
give the asymmetry required. The values required would not lead to 
convenient waveform generation; however, it is not our purpose here to 
propose a practical design. 

Figure 5.4 shows the potential profiles obtained applying the program 
to a silicon substrate. TSI, QUN, SIGMA, and YBAR are chosen to 
obtain an N-layer 1 microns thick with a doping of 3 x 10^^cm“^. The 
P-substrate is doped with lO^^cm”^. The resultant profile in A would 
not allow undirectional charge transfer. Some items of interest for 
this curve are as follows. The residual dropped from 4.52V to ImV in 
240 iterations requiring 149 seconds CPU time. The depletion depth 
occurred at y = -5.25y and h = 0.48p giving a maximum roundoff error 
of approximately 10%. The field contains 4047 grid points. 

The curve labeled "B" shows sufficient asymetry to make charge 
transfer to the right possible. The maximum potential occurs at a 
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distance from 0 to 0.4u below the Si-SiO^ interface. Unfortunately, 
under the pad where charge would be transferred to, the maximum occurs 
on the surface. However, these results are given to illustrate the 
program capacilities and not to suggest a practical design. 


5.3 Use of Program 

The program is approximately the same size and takes approximately 
the same running time as the 3-electrode program discussed in Section 4. 
Convergence is obtained within 240 iterations to within 3% or better with 
the residuals typically dropping from approximately 5 volts to lOmV. Such 
a run requires 158 seconds with 4047 grid points in the field. 

The dimensions of the arrays must be checked for compatibility with 
the desired number of grid points. The dimensions in the listed program 
allow 100 X 100 points. The listing is given in Appendix C. 

Input data cards are in the sequence given below: 


data format 

1. LP1,LG1,LP2,LG2,1TYPE, MASCO 6110 

lTYPEs=0, sapphire substrate; ITYPE = 1, silicon. 

MAXU=1 for printout of u , otherwise MAXU=0. 

' max 

2. JA,JB,JG 3110 

JA: lower line printed out; JB: upperline; 

JC: number of lines skipped. 


3. L1,L2,L3,L4 


4110 


4. L5,L6,L7,L8 


4110 


5. TOX1,TEL1,TOX2,TEL2 4F10.3 

6. T0X3,TSI,W 4F10.3 


7. QSS1,QSS2,QUN,CSUB 4E10.3 

QS^l: Si-Si02 interface surface state density, cm“^. 

QSS2; Si-Sapphire Interface surface state density, cm“^. 
QUN: Uniform implanted density, cm“^. 
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DATA FORMAT 

8. VGl, VG2,VG3,VG4, ITER, LAP 4P10.3,ZI10 

(see Section 4 for more Information.) 

9. OMl,OM2,OM3,OM4, YEAR, SIGMA 6F10.3 

Try omega value between 1.5 and 1.8 

Subroutine DISFAC uses YEAR and SIGMA as 
the usual parameters of a gaussian function in 
coordinate variable y to specify implanted doping. 

10. - Repeat card (8) for as many values as desired. 

(Last) Blank card terminate program. 


References 


1. A. B. -M. Elliot and J. C. Anderson, Solid State Electronics, 
Vol. 15, pp. 531-545, 1972. 
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6. CONCLUSION 

We have described three programs which can be used to obtain an 
electrostatic analysis of CCD structures. These programs cover a broad 
variety of structures which may be of interest. When the computed pro- 
files are such that uni-directional charge transfer is possible, then 
such profiles may be considered to be approximately those which would 
exist in a structure operating in an ideal mode with a small amount of 
"signal" charge in any given cell. The profiles also serve as a base 
for estimating the amount of charge corresponding to a "full well". 

Other work which we have done on CCD’s during the past year includes 
dynamic analysis of charge transport. Our work has been confined to one- 
dimensional models, and we have studied the so-called "flux-corrected" 
method for solving the transport equations. We are not ready to report 
on this work at this time. The emphasis which has been given to buried 
layer devices and the observations which we have made from electrostatic 
analyses makes it clear to us that accurate models of the transport 
equation must also be two-dimensional. If one does not use two-dimensional 
models , then it seems that the techniques used by Kosonoeky and Carnes 
and Amelio (references in Section 1) are adequate for estimates of the 
transit time. It should be noted that a two-dimensional analysis will 
require at least twice as much data storage as the electrostatic 
analysis. Computing time will be at least doubled and probably increased 
well beyond that. However, we intend to pursue this problem further. 
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APPENDICES 

The following listings are of programs written in FORTRAN V. 

All programs are approximately the same size. Program 2, Appen- 
dix B, requires 20,000 words or 80,000 BYTES of storage. The field 
points require approximately 11,000 words. Without changing the arrays, 
all programs handle a grid of 100 x 100 points. (Appendix B program 
100 X 110) . Runs will require up to 240 iterations per solution with 
about 4 minutes on a UNIVAC 1106 , or equivalent , for one solution 
using 10,000 points. Reduction of the number of field points will 
reduce the run time proportionally. 

To use, supply run control cards, with appropriate control 
cards also separating the main and subprograms , and prepare data cards 
describing structure according -to formats given. , Suggested trial 
values of relaxation parameter are between 1,5 and 1.8 with 1.7 a 
good value to start with. 
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APPENDIX A 

THREE-ELECTRODE PROGRAM 


See section 3, Figure 3.1 for diagram of structure. 
Data cards : 


DATA FORMAT 

1. JA,JB,JC 3110 

JA: Lowest line of potential array printed. 

JB: Highest line of potential array printed. 

JC: Spacing between lines. 

Program recompute JA and JB to obtain printout 
of line J2, the Si-S102 interface. 

2. TOX,TEL,W,CDEBY,CDEPL 5F10.3 

Default values of CDEBY and CDEPL are 2. 

See text for discussion of optimizing these values. 

Program will run well on default values. 

3. QSS,CSUB,QG,QP 4E10.3 

See listing and text for further information. 

4. V1,V2,V3,ITER,LAP 3F10.3,2I10 

Try first runs with ITER = 240, LAP = 40. 

This will give approximately 50 pages of output 

data including listing. After inspection for convergence, 

both ITER and LAP may be reduced. 

5. L1,L3,L4,L5 4110 

See Figure 3.1 text. For example, 

LI = Jl-1 , L2 = J2-J1, etc. 

6. LP,LG,LDOPE,LSPAC 4110 

See listing and text for further information 

7. Repeat card 5 for as many solutions as desired. 

8. Last data card is blank to stop computation, 



»RODUCIBILITy OF THE 
O^amAL PAGE IS POOR 
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X.DECKf »MINE 
B 6 9_ 01/28"! i; 47 ;17 


..—Ofll) oauu -PCD^i£EOOl7UjLJljs/j+fy5o 

000 OASG»T TEMP»F2 

OfliJ RFLTfL.lB TFM P.PCH 

Obb C *****st!*********5*!!t: + + ***:*!***j(;:t.*^** + + :k + **;4j* + ***!(c*4r******+****:((** 

- _0P0 C ELE.C.TRQSIATI_C_ ANALYSIS OF CCD STRUCTURE 

000 c dimension arrays » define common variables 

000 c H is the ePXD_SPACljig_NQfmAL Tn THF iMTrRFACF 

000 c w » * ♦ Parallel with the interface 

000 c u IS the electrostatic potfntial wrt the substrate bulk 

OOp C DU IS the residual 

000 C ***** ♦****♦*** ♦ **♦»♦* *♦***»* * ;****.** * ************************* 

000 COMMON U(lOOflob) f dU( lOO ) »H C lOO ) »F0»L2, ALPHA » IMAX »JMAX » EA rEB » 

QOO *Ajr.B.f..C.t^Wjg..V_T t Ji.DUM AX_ 

00,0 DIMENSION QS<100) 

QOO .VT=0.0259 

000 c*****************?*"************************************************* 

OOP G READ LINE _mi NTQUT, CONTR OL Pa RA MFTFRS 

000 C**=t‘***^*******************’<‘*******4;************************************ 

_000 REM15jJJ_.ua » JH ». uc . : 

000 1 format ( 31 lb) 

■ 000 C ♦***************»************,|.******»*A****»***************** 

000 C READ DIMENSIONS OF REGION FOr ANALYSIS# MICROnS 

000 C ********** * ** ****** *****!»^t* ^c*:fc****** ******* ♦ ♦ ****♦*♦***♦♦♦*** 

000 c *****************************^******************************* 

boo C T0X= OXIDE thicknfss 

000 C TEL= ELECTRODE THICKNESS 

000 READ (Sr 2) TOX >TEL # W # CDEBY > CDEPl . . . 

000 2 FORMAT (5F10. 3) 

QOp C QSS= SURF ACE STATE CHAR GE > C M-2 ^ 

000 e CSUB= SUBSTRATE DOPING rCM-3 

boo C QG= GAP BOPI NGrC M-grPO SITlVE TF CHAR GE POSITIV E ^ 

000 C QP= UNDERPAD DOPING' * ♦ »»• 

000 C *****************************ifc***********»*»***********»***** 

000 REA0(5r4)QSSrCSUBrQGrQp 

000 4 FORMAT(4E10,3) 

poo c ’SIGN' IS used To PRECLUd"e s'paCE CHARGE INSTABILITY IN RELAXATION 

000 C IN SUBSTRATE . IF *U t HAS WRONG STGNr »U« IS SET TO ZERO. 

000 SIGN=-1*0 

000 IF(CSUB.GT.O.O) SIGN = l.n ^ , 

000 C ♦♦*»♦♦*♦♦********♦**♦***♦*♦***************♦********♦********* 

000 C READ ELE CTRODE VOLTAGESrNUMRER oF ITERA TI ONS# A MD LO OP PRINTOUT 

000 C ************************************************************* 

_ pop READ (5»6) V1#V2 # V 3 rJTERj LAE. 

000 6 FORMAT(3Fl6,3#2il6) 

QQO C =♦;***♦*********»***»♦**♦*♦♦****♦*♦♦*♦»****»★♦*»*♦***♦*»*»*♦**♦ 

poo C COMPUTE CHARGE DENSITY PARAMETERS 

boo C **** * *** ***** *** * *^ * ***j^ ********************************** 

000 (3=1,6E-19 ■ 

OOP E0=ft.R5E"14 • 

OOQ VSS=Q*QSS/E0*l,E-4 

OOP VSUR=«Q*CSUB/EQ*1 .F"8 

000 C0N=2,*11.7*E0/(Q*CSUB) 

.. 000 C * ***** ******* *******»**** ****^**i){*****^*i4e****^(^*:jc*»4;*»****ite»i)t 

000 C ESTIMATE DEPLETION DEPTH ' "" 
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umj — 
odo 

c 

c 


ESTIMATE SURFACE pITTEnTIAL 

000 

000 



ABV1=ABS(V1> 

ABV2=ABS(V2) 

000 

000 



ABV3=AB5(V3) 

VG=V1 

000 



IF(ABV2.GT.ABV1) VG=V2 

000 



rFU8V3.GT.ABV2> VG=V3 

000 



C0=3. 9*8. 85E-1A/ ( 1 ,E“^*TOX ) 

000 



VA=VG+QSS*Q/CO 

000 



VB=Q*CSUB/C0*11 .7*eo/co 

ocio 



PSI= VA+VB“SQRT(2.*VA*VB+VB**2) 

000 



XDS=2 , *1 1 . 7*E0/Q*PS I /CSUB 

ooo 



XD=SQRT(ABS(XDS) )*1.E4 

000 

c 


«4;4:ic4c«3|:*^4:*4:**«*3(c****3(c*4:«*4c:4c*^4:«4:4:****4c************^********* 

000 

c 


FIND THE EXTRINSIC DEbYE LENGTH 

000 

c 


4t4c4c:^3te:|c:(c4E«4c*4:**4c**************:^***4:****#********************** 

000 



DL=\/T*11.7*E0/Q/CSUB 

000 



EDL=SQR-T(DL)*1.E4 

000 

c 



000 

c 


WRITE DEPLETION DEPTH AND DEBYE LENGTH 

000 

c 


4i)|c4:«*^:4t*4t****it:4:***«***«*4:**:4‘*4;*4:3(c3t:****4c**4:***3tc:tc******4c*****4:* 

000 



WRITE{6»300> XD^EDL 

000 


300 

FORMAT(’l»rlOX» 'DEPLETION DEPTH =♦ »F10,3»5X» 'DEBYE LENGTH=» » FlO . 3) 

“OOC5 

c 


*4c*4:*********«»***4:**»:»:**«***4;***********4c******4c**«*****4:*4c* 

000 

c 



odd 

c 


determine normal control indices 

000 

c 


AND SET THE Y-GRID SPACING 

000 

c 


Ll» REGION OF UNIFORM SPACING IN SEMICONDUCTOR 

000 

c 


L2»» ♦ EXPONENTIAL • » ♦ 

IRTO 

c 


L3» » » UNIFoF^M ♦ T“ OXIDE 

000 

c 


L4r » » ♦ ' through electrodes 

000 

c 


L5» • ' EXPONENTIAL * * IN AIRSPACE 

000 

c 


♦ ♦****************************:*;»**!':**********+************1t** 

000 



REAO(5»302) Li?L3rL4»L5 

000 


302 

FORMAT(4I10) 

000 



J1=L1+1 

000 



K1=L1 

000 



IF(CDEPL.LT.i.d) CDEPL=2 .o 

000 



DO 8 J=lfKl 

000 



H(J)=CDEPL*XD/Ll 

000 


8 

CONTINUE 

000 



H0=3.*T0X/L3 

000 



IF(H0,6T.EDL) HO=EdL 

000 



IF(CDEBY.LT.2.) CDE8Y=2. 

000 



Y0=C0E8t*Hd 

000 



FO=YO/HO 

odo 



CALL YSPAC 

000 



J2=J1+L2 

000 



K2=J2-1 

ODO 



DO 10 JJ=lrL2 

boo 



vJ=J2-JJ 

00.0 



H(J)= HO*ALPHA**jj 

000 


10 

CONTINUE 

000 



J3=J2+L3 

000 



K3=J3-1 


000 ■ ■■ DO 12‘j=J2rk3 
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iia 


CONTINUE 

_U5=Jit±L5. 

K5=J5-1 


H(J) = (TEL/L4>a(2«**JJ1 

16 CONTINUE 

♦ »♦*♦** « ♦ *»* * ♦ * * »-f* *** ** * * ****** * *** ************* 

PRINT OUT THE Y-COoRDINATE VAlUeS AND THECONTROL INDICES 

******** ***** * * ******* * * jM: ** *jfe* ****** * ******** ** ** * * **** * **** 

WRITE C6» 17) 


Y=0. 

DO 12S J=J2 >K5 

Y=Y+H( J) 

WRITE(6»1 27) Y ^ 

125 CONTINUE 



DO 129 LL=1»K2 

J^2d-L„ _ _ 

Y=Y-H(J) 

WR ITE (6» 127) Y 

129 CONTINUE 

127 F0RMAT(F2Q.6) ^ 

WRITE(6»131) Jl» J2f J3» J4» J5 

131 FORMAT(// »10 X> *J"VaLUES ARE! * , 5Tl n//> __ 

c compute ja»jb values to obtain interface printout 

MT=(U 2~JA )/JC _ _ 

JA=J2-MT*JC 

IF(jBtGT«»J5 ) vJB~U 5 , 

:*:**^c********** + ***!(t****:**!** + *^*******:***3(r***** + *******c*****:*!* 

C DETERMINE LATERAL CONTROL iNplCES 

C LP AND LG'aRE THE NUMBER OF STEPS ACROSS THE PAD AND GAP. 

^ T HEN CO MPU TE THE SURFACE CilA.RGE. DISTRIBUTION 

C **♦** + ♦*♦ *.******** + *******:***:4t******************************* 

READ(5>18) LP»LGrLpOPE.LSPAC 

18 FORMAT (41 10) 

I^li+LP' ■ 

I3=I2+LG _ _ _ 

i4=l3+LP 

■I5=I4*LG : 

I6=I5+LP 

I7=I6+L P __ _ ■ 

lA=(I2+I 3)/2 

IR =(i 4 -K 5 )/2 REt'KUD U Uib ■ I .ITY QF~ TH H’. 

ID=lC-<-LDOPE QRIGIHAL PAGE IS POPE 

IE=I 3 +LSPAC 

If =JE+LD 0 PE _ _ 

IG= 15 +LSP'aC 
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— crono— 
000 

■■■■oO"d"' 

000 

■ - 

- 

IH=r5+CDT5PET" 

DO 333 1=2 » 17 

■ QS ( i ) =QSS+QG* (F(lfIl»l)+F(I2»l3tI)+F(I4»'T5»I)+F(l6»I7»I)) 
♦+QP*(F(IC»ID»I)+F(IE»IF»I)+F(I6»IH»I) ) 

000 

000 


■ ‘q'S (11) =QS (ID rQG/2 . 

333 CONTINUE 

000 

000 



QS(I2)=QS(I2)-QG/2. 

QS(I3)=QS(I3)-QG/2. 

oob 

000 



Q5(I4)=Q5(I4)-QG/2, 

QS(r5)=QS(I5)-QG/2. 

000 

000 



QS ( 1 6 ) =QS a 6 ) -q'G/2 , 
QS(IC)=QS(IC)-QP/2, 




QS(ID)=QS(ID)-QP/2, 

QS(IE)=QS(IE)-QP/2, 




QS ( IF ) =QS ( IF ) -QP/2 , 
QS(IG)=QS(IG)-QP/2. 


c 


QS(IH)=QS(IH)-QP/2, 

*****:f:****s|:***********!»:*itc****************************#*****>t:* 


c 

c 


SET THE ELECTRODE POTENTIALS 

*♦*♦*♦*♦♦****♦*♦**♦♦♦*****♦***♦**♦****♦♦*♦*♦♦♦♦♦**♦*♦♦♦♦♦**** 




DO 28 J=J3rU4 
DO 26 I=I1»I6 



26 

U d » J ) =V 1*F ( 1 1 » 1 2 » I ) +V2*F ( 1 3 » 14 » I ) +V3*F ( 1 5 d6 » I ) 
CONTINUE 


c 

28 

CONTINUE 

000 

000 

c 

c 

SET VALUES BETWEEN ELECTRODES 0^t AIR-OXIDE INTREFACE 

♦ ♦♦♦♦**♦♦*♦*♦♦♦**♦★*♦***♦*♦♦♦+♦*♦*****!♦:♦♦♦*♦***********♦♦♦★♦♦ 

000 

000 



J=J3 

DO 240 I=2rll 


240 

U(If J)= U(lrJ)+( I-l)*<U(llf J)-U(lr J) )/(Il-D 

CONTINUE 




"DO 242 1 = 12 D3 

Udf J)=U(I2»J)+(I“l2)*(U(I3f J)-U(I2»J) )/(l3-I2) 



242 CONTINUE 

DO 244 I=I4»I5 


244 

U(I»J)=U(I4»J)+(I-I4J*(U(I5»U)-U(I4» J) )/(I5-I4) 
CONTINUE 




DO 246 I=I6»I7 

Udr J)=U(I6»J) + (I-I6)*(U(I7» J)-U(l6r J) )/(l7-l6) 


C 

246 

CONTINUE 


c 

c 


INITIALIZE POTENTIAL DISTRIBUTION IN THE SEMICONDUCTOR AND 
OXIDE USING THE ONE DIMENSIONAL DEPLETION LAYER ESTIMATE. 

000 

000 

c 

c 


VB DEFINED ABOVE 

♦ **♦*;*{***♦***♦:+:♦****♦**♦**#*♦*****♦*♦♦*****♦*♦**♦♦**★♦♦♦♦♦*♦* 

000 

000 



DO 406 1=2 » 17 
J=J2 

000 

000 



VS=U(I»J3) 

VAA=Q*QS(I)/CO • - 




VA=VS+VAA 

IFd.Gt.l.AND.I.LT.Il) GO TO 39.4. 




IF(I.GT.I2.AN0.I.LT*I3) GO TO 398 
IF(i.GT.I4.AND.I.LT.I5) GO TO 398. 




IF(I.GT.I6.AND,I.LT.I7) GO TO 398 
VlNT=VA+V8“SQRt ( 2 , ♦VA*VB+VR**2) 




GO to 399 
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■■ 398 
399 

V INT=Q* IQS(I)**2)/(2.*11, 7*E0*CSU8 > 
.tONTJl^UEL 

OOQ 

000 


XDS=CON*VINT 
.Xb=SQRT (ABS(XDS) ) 

000 

000 

400 

Y=0, 

J=J-1 

000 


Y=Y+H{J) 

000 


_ YN=1».E-4*Y/XD 



IFCYN.GE.I.) GO TO 402 

U(I»J)=\/INT*{1 .-*Yn^**2 


402 

GO TO 400 
J=K2 

000 


Y=0. 

_-Q00„_ 

404 

J=J+1 

000 

-j3m .. 


Y=Y+H(J) 

YN=Y/TOX 

000 

000 


IFCYN.GE.l.) GO To 406 
U(I,J)=VINT+{VS-VImT)*YN 

000 

000 

406 

GO TO 404 

__£:iMTJ.NUE _ 

000 

000 

C WRITE OUT THE STRUCTURE PARAMfTFRS 

000 

000 

C^*********#*** + **************))c3(c#:f:*^)j:**+*+ + *3(.:4(*:j:***:j!3jt****5t:*********^:f.**** 

WRITE(6>345). T0XyTEL»WrVl»Vg»v3»LPyL6rLD0PE>LSPAC»Qy;g;»QG»QP»rg;UR 

ooo 

000 

345 

FORMAT(/»10X» ♦TOX»TELrW=:’ y3Fl0.3/10X»*VlfV2rV3=» »3F10.3/ 

*10 X y ♦_LP » Lj3.y LDOP.Ey LSPAC= »f4IlO/lnX»QSSyOGyQPy CSUR= » y 4E1 5 . 3/ ) 

000 
000 
000 
00 0 

C 

._C 

C 

♦ »******#********************^:*;****************<t*****)»:******* 

_ R£ A D_ RELAXATION PARAMETER VALUER 
READ{5y408) 0Mly0Mgy0M3y0M4y0M5 

000 

ood 

408 

FORMAT(5F10.3) 

-_kt=l 

000 

000 


JL0=Jl/2 

JHI=J4+1 

000 

000 

C 

C 

♦ ♦♦♦♦:«!+:******** + ******4!******** + ****+*^*********^****»******* 
START THE RELAXATION PROCEDURE STARTING AT THF LOWEST ROW 

000 

doo 

c 

c 

AND WORKING ACROSS TO THE RIGHT TO THE TOP ROW 

♦ ♦♦♦♦*#*********#**** + ***********^c*************:*:*^***:4c**:^^**^ 

000 


LOOP=l 

000 

800 

CONTINUE 

000 


DUMAX=0. 

000 


DO 40 J=JLOyK? 

000 

000 


EA=11.7 

EB=11.7 

000 


A=EA*W/H( J-1) 

000 


B=(EA*H( J-1)+EB*H( j) )/(2,+W) 

000 

000 


G=EB*W/H(U) 

C=A+2., ♦8+G _^^T.7TnTO7T.TTY OF THE 

000 

000 


OMEGArOMl PAGE IS POOR 

DO 38 1=2 y 17 OBIGIN^ 

000 

000 

30 

IF(UnyJ+l) .EQ.0.0) go t’6 33 
-J'JE.STSABS ( U ( I y J-1 ) /VT ) 

000 

000 


IF(TEST.GT.4> GO TO 31 
.ITEST=1 

000 

000 


QI J=VSUB*AVERO ( ITEST) 
..GQ_T0 32 


000 31 QIJ=VSUbVw*'(H( J)+H( J-1) )/2”o 
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00 0~ 
000 

32“ 

UOLD=Ua»J) 

UTIL=(A*U{I»J-D+B*U(l-l»J)+B*U(I+lf J)+6*U(I»J+1)+0IJ)/C 


000 

000 


U ( I » J ) = { 1 . -OMEGA ) *uOLD+OMEGA*UTIL 
US=5IGN*U(I»J) 


000 

000 


IF(US.LT.d.O) Ud»j)=0,0 
nu{I)=U(I»J)-UOLD 


000 

000 

33 

GO TO 36 

IF(U(I+1» J) .EQ.0.0) GO TO 34 


000 

000 

34 

GO TO 30 

IF(U(I-i»J) .EQ.0.0) GO TO 38 


000 

000 

36 

38 

CALL MAXOF 
CONTINUE 


000 

ood 

40 

c 

CONTINUE 


000 

000 

c 

c 

RELAX the SILICON, silicon-dioxide INTERFACE 

♦ ♦♦♦♦♦♦♦♦★i^t****^**************;^:***************************^^** 


000 

000 


J=J2 

EA=11.7 


000 

000 


EB=3,9 

A=EA*W/H(J-1) 


000 

000 


B=(EA*H<J-1)+EB*H(J) )/(2.*W) 
G=EB*W/H(J) 


000 

000 


C=A+2.*8+G 

0MEGA=0M2 


000 

000 


DO 48 I=2fl7 
GSA=(QS(I+1)+QS(I-1> )/2. 


000 

000 


fESt=A8S ( U ( I f J-1 ) /Vt ) 
IF(TEST.GT*4) GO TO ‘*2 


000 

000 


ITEST=3 

QI J=VSU8*AVER0 { ITEST) +Q*QSA/E0*1 *E-4 


— mr 
000 

42 

GO TO 44 

QIJ=VSU8*H(J-1)»W/2.+G*QSA/E0:<c1,e-4*W 


ooo 

000 

44 

UOLD=U(I»J) 

UT I L= ( A*U ( I » J-i ) +B*U ( I -1 r J ) +B*U ( I+l » J ) +G*U ( I f J+1 ) +QI J ) /C 


■ffOlT" 

000 


u a , j r^r. -omegai *uol'd+omega*uti l 

US=SIGN*U(I»J) 


oUo~ 

000 


IF(US.LT.O.d) U(l,jr:^.0 
DU{I)=U(I,J)-UOLD 


000 

000 

48 

CALL MAXOF 
CONTINUE 


000 

000 


DMAX=0.0 

J=J2 


000 

OOP 


DO 696 1=2,17 
IF(ABS(DU(I) ) .GT.DmAX)IM=I 


000 

000 

696 

IF(ABS(DU(I) ) ,Gt.DMAX)DMAX=ABs(DU(I) ) 
WRITE(6,697 )DMaX,Im 


000 

000 

697 

C 

FORMAT ( 2X, •DU=» »F5.3,5x,»Ir» »i3) 


000 

000 

c 

c 

RELAX THE OXIDE REGION 

- 

000 

000 


JJ=J2+1 
DO 60 J=JJ,K3 


000 

000 


EA=3,9 

EB=3,9 


000 

000 


A=EA*W/H(J-1) 

b={ea*hcj-i)+eb+h(j) )/(2.*W) 


000 


G=EB*W/H(J) 





67 


000 

000 



C=A+2.*B+G 

0mEGA=0M3 

000 

000 



DO 58 I=2»I7 

Qij_=a.o_. . _ _ _ _ 

000 

000 




000 

000 



U ( I » J) = ( 1 . -OMEGA ) *UOLD+OMEGA*UTIL 
DU(T)=U{Ir J)-UOLD 

000 

000 


58 

CALL MAXOF 
CONTINUE 

000 

000 

_JL 

60 

CONTINUE 

000 

000 

c 

q 


RELAX THE OXIdE-AIR-ELECTRODE INTERFACE 

♦** + *»***^55t:*)*:*)(c***+***)tc******* + ***!*:******^:*:J:*****:(t*****:ir*****: 

000 

000 



J=J3 

EA=3.9 

000 

000 



EB=1.0 

A=EA*W/H(J-1) , , 

000 

000 



B=(EA*H(U-1)+EB*H(J) )/(2..*W) 
G=ER*W/H( J) 

000 

000 



C=A+2,*B+G 
DO 68 I=2fI7 

000 

odo 



IF(I.GE.I1.AND.I.LE*I2) go TO 68 
IF( I.GE.I3.AND.I.LF.I4) SO TO 6R 

000 

000 



I F ( I « GE • I 5 • AND « 1 • LE • ) GO TO 68 
QIJ=0, 

OOO 

000 



0MEGA=0M4 

UOLD=U(I#J) 

000 

000 



UTIL=(A*U(I»J-1)+B*U(I-1» J)+8*U(I+1» J)+G*U<I»J+l>+QId)/C 
U { I mJ) = { 1 . -OMEGA ) *U0LD+0MEGA*UTTL 

000 

000 



DU(I)=U<I»J)-UOLD 
CALL MAXOF 

odo 

000 

1._C. 

68 

CONTINUE 

000 

odo 

c 

c 


RELAX THE ELECTRODE-AIRSPaCE REGION 

000 

000 



JJ=J341 
DO 80 J=JJfJ4 

OOD 

obo 



EA=1. 

EB=1. 

000 

000 



A=EA*W/H( J-1) 

.B=(EA*H(J-1)+FB*H(j) )/(2.*w) ... 

000 

000 



G=EB*W/H(J) 

C^A+a.^t'B+G 

000 

.000 



OMEGA=OM5 ^EPIiODUCrBILlTY OP 

DO 78 1=2 » 17 ORIGINS T. PAGE JO 

000 

000 



IF(I.GE.I3.ANd.I,.LE»I‘<-) so to 78 ^ 

IF{I.GE.I5.AND.I'.LF*I6) GO TO 7.8 • ' 

000 

000 



IF(I.GE.I1.AND.I.LE»|2) go to 78 
.010=0. 

000 

000 



U0LD=U(I»0> 

UT_IL=(A*U(I»J-l)+B*U{r-l»J)+B*U(I+lr J)+G*ua»J+l>+OIJ)/C 

000 

000 



U < I r J ) = ( 1 . -OMEGA ) ♦UOLD+OMEGA*UTIL 
DU(I)=U(I» J)-UOLD 

000 

000 
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CALL MAXOF 
CONTINUE 

000 


80 

CONTINUE 
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Tjinr 

^oo_ 

000 

_PM, 

000 

000 


c 

c 


“♦»*♦***♦♦♦♦*♦♦♦***♦♦♦*♦♦***♦**********♦*♦**************♦♦***♦ 

_ REL^_ THE air- space REGI ON A BOVE T HE ELECT RODES 

**♦♦**♦♦*♦♦♦♦!♦;♦*♦*♦♦**♦♦♦*♦♦*♦***♦**♦****♦*♦♦*♦♦*♦♦*♦♦♦**♦**♦ 

_ v)J=J4+l 

DO 90 J=JJ»JHI 

EA=1* 


000 

0.00 

000 

000 


EB=1. 

A=EA*W/H ( J-1) 

B=(EA*H(J-1)+EB*H(J) >/(2.*W) 
G=EB* W/H(J) 


000 

000 


C=A+2.*B+G 

0MEGA=0M5 


000 

000 


000 

000 


82 


DO 88 I=2»I7 

QIJ=0. 

IF<U(I»J-1) .EQ.O, ) GO TO 84 
CONTINUE 


000 

000 

000 

000 


UOLD=U<ir J) 

UTIL=(A*U(I>J-l)-*-B*U<I“l»J)+B*U(H-l»J)+G*U<I»J-H)-*-QIJ)/C 
U ( I » J ) = ( 1 • -OMEGA ) *U0LD+0MEGA*UTI L 

DU(I)=U(Ir J)-UOLD 


000 ■ CALL MAXOF 

OOP GO TO 88 __ 

000 ■ 8.4 IF(U(I+lr J) .EG.O. ) GO TO 86 

000 GO TO 82 

"WO EE IF(U(I-l»Ji .EQ.O. ) GO TO 88 

000 GO TO 82 

000 88 CONTINUE 

boo 90 CONTI_NUE 

000 ■■■ ■ ULOrJllb-l 

boo IF(JLO.LT.2) JL0=2 ^ 


000 

000 

iroff" 

000 

000 

000 


GOO JHl=JHl+l 

000 

c ' ’ :^iiic*;ir**;(c*i*:**** + ***********^ + **^ + **^e*******)t:***^********* + ***** 

c OUTPUT control AND ITERATION CONTROL FOLLOWS _ 

rEPLA^D"BY X CONDlTlbNAC 

C CONTROL ON either QU OR DUMAX, 

C ♦**»♦*#♦♦**♦*♦*♦♦*♦*♦«*♦♦*#*♦*:♦£♦♦*****♦*♦*****♦*♦*♦**♦****♦*♦ 

LOOP=L OOP-H . 

00 0 "■ LWRIT=LOOP/LAP 

000 IF(LWRIT.EQ.KT) GO TO 881 ^ 

000 207 ' CONTINUE ’ 

000 ^ IF(L 00 P»LT.ITER) go to 800 ^ 

00 0 C**:t:**:t:*******#************Jf;*******^*********!tt*******+**i»t***********+*** 

000 c The follow.ing read may be used for finding new potentials ^ 

boo C ■ after applying small "perturbations to THE electrode VOLTAGES 

000 r ************************************************** ********************* 
‘oob RE^( 5 * 6 ) Vlf‘V 2 »y 3 ,ITERrLAP 


u w 
000 
000 


i r ' A » urv ♦ UV« • ^ 1 

DO 883 J=J3»J4 
DO 882 I=I1»I6 




U(Ir J)=Vl*F(li»i2»l)+V2*F(I3»l4rI)+\/3*F(I5>I6»I) 
♦+(l.-F(Il»I2»I)-F(i3»l4»I)-FCi5»I6»I>)*U(I»J) 


000 

■ 882 

CONTINUE 


000 

883 

CONTINUE . 



000 

OpO 

bob 


L00P=1 

KT=1 


GO TO 800 
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000 

000 

881 

CONTINUE 

WRITE(6»212> nUMAXtIMAXf UMAX 

000 

000 

212 FORMAT ( 30X » *DUMAX=» f FlO .3» 5X » » IMAX=» » I4» 5X» ♦ JMAX=* » 14) 
no ?08 J=JA»JR*JC 

000 

noo 

HR 

WRITE(6»208) J 
FORMAT(10X»'J=»»I3). 

boo 

000 


WRlTE(6f206) (U(IfJ)» 1=1, lA) 
WRITE(6,206) (U(I,vJ)» I=TA.TR) 

000 
000 . 

__ao_6_ 

WRITE(6,206) I = IB,I7> 

000 

000 

205 

CONTINUE 

_KT=KT+1 

000 

000 


GO TO 207 
END 

000 

000 


FUNCTION F(II,JJ»KK) 
F=0.0 

poo 

000 


IF{KK.6E.I1.AND.KK,LE,JJ) f=i.o 
RETURN 

000 

boo 


end 

SURROI'JTINF tSPAC 

000 

000 


COMMON U(lOOrlOO) , dU ( iOO > , H ( iQO ) ,FO»L2 , ALPHA, lMAX»JMAXfEA,EB, 
*A,R»C.G,W,VT,I,J,m IMAX 

000 

000 


ALPHA=1,25 

L2=AL0G ( 0 . 25*F0+1 . 0 ) /AlOG ( ALPHA ) 

000 

000 

2 

Fl= { ALPHA**L2-1 . ) / { ALPHA-1 . ) 
TRY=ABS(F1/F0) 

000 

OPO 


IFtTRY. LT, 1.005 . and. try. GT. 0.995) GO TO 4 
F=F1-F0 

000 

obo 


LL=L2-1 

F3=L2*(ALPHA**LL)/(ALPhA-1, ) 

000 

000 


F4=(ALPHA**L2-1. )/( ( ALpHA-1 . ) **?) 
FP=F3-F4 

000 

000 


ALPHA=ALPHA-F/FP 
GO TO 2 

000 

oob 

4 

CONTINUE 

RETURN 

000 

opo 


END 

SUBROUTINE AVEROdjEST) 

000 

000 


COMMON U<100,100) »OU( lOO) ,HdOO) »FO, L2, ALPHA, lMAX»JMAXfEA,EBi 
♦a,b,c,g»w»vt,i,j,dumax 

000 

000 


5=0.0 

Ul=ABS(Ud,J)+UdrJ-l) ) 

000 

000 


IF(Ul.GT.O.O) S=0.5*H(J-1) 
IF(ITEST,6T.2) 60 TO 4 

odb 

000 


U2=ABS(Ud»J)+Ud,J+l) ) 
IF{U2.GT.0.0) S=S+0.5*H(J) 

000 

000 

4 

aver’o=s*w 

RETURN 

000 

000 


END 

SUBROUTINE MAXOF 

000 
op a 


COMMON U(100,100) ,dU( 100) ,H(l00 ) »FO,L2 , ALPHA, ImAX»JMAX»EA,EB, 
*ApB , C , 6f W , VT r I , J , OUMAX 

bob 

000 


D0LD=ABS(DUMAX) 

DNEW=ABS(DU(D) 

000 

000 


IF(DNEW.GT.DOLD) GO To 2 

_..GO TO 4. . _ _ 

000 

2 

DUMAX=DUd) 




"GOO"” 

000 


IMAX=I 

JMAX-J 





000 

4 

RETURN 





000 


END 





000 

qxqt 





- 

000 


5 40 


5 



poo 

0.5 

0.8 

2.0 




000 

2.00E+11 5.00E+14 




*• 

000 

0.0 

10.0 

16,0 


240 120 


000 

• 

25 5 


5 

10 


000 


8 4 





000 

1 

.8 1.8 

1.5 


1.5 1.5 



0 


(SEND 

flPCHfS TEMP.PCH 
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APPENDIX B 

TWO-ELECTRODE PROGRAM 

See Section 4, Figure 4.1 for diagram of structure. See end 
of Section 4 for description of data cards. End of listing gives 
example . 
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" Me® IS POOR 


X*DtCK» rWTTJE 

Ql/a8"li:49t52 


000 QRUN CCD 2 »EE00l7 Ur JIMr r/50 

000 QASGrT TEMP?F2 

000 QELtrLiR TEMP.PCH 

000 C;*:#***4:** + ******************»****!t:*!(e*********************************** 

dob C ' ' PROGRAM MUMBER 2******PERIQDlC BOU NDARY CONDITIONS************* 

' boo c**’*^*********** TWO OXIDE thicknesses ARE ALLOWABLE********************* 

QQQ Q A A *:fc* A***** **=*********** ***************** * ******************* 

000 c electrostatic ANALYSIS OF CCD STRUCTURE 

000 G DIMENSION ARRAYS » DEFINE COMMqN VARIABLES 

000 e H IS THE GRID SPACING NORMAL TO THE INTERFACE 

000 c w » ’ ' Pa rallel with the interface ^ 

000 C U IS the ELECTROSTATIC POTENTIAL WRT THE SUBSTRATE BULK 

000 C DU IS THE RESIDUA! 

boo C ****#******|*:*i*********************************************** 

OOQ COMMON U(100»110) yPUtlO O) yHdOOi >F0»L2 » ALPHA t 

000 ♦EArEBf A»BrC»G»W» VT rSIGN»DROW»DUMAX» ABMAX» 

000 * I f J » lM AX»JMAXylM»l5 

"" 000 ■ DIMENSION QS ( 100 ) » CHI0n( ll ) » ARRflY ( 100) 

000 SIGMA=0.1 

000 YBAR=-0.2 

boo VT=0,Q259 ^ 

Wo READ(5»1) JA» JBrJC»MAXU»lRES 

000 1 FORMAT(5I10) 

blTD C ♦♦♦iicljt^****************************************************** 

000 C read DIMENSIONS OF REGION FOr aNALYSI S> MICROnS 

— oW C ************************************************************* 

000 C T0X1= THIN OXIDE THICKNESS. 

5W C TOX2= THICK OXIDE THICKNESS 

000 C TEL= ELECTRODE THICKNESS 

Obb"" “ ” READT5»2) T6Xl»TOX2»TETLrW»HFAC»HbEBY 

000 2 FORMAT (6F10. 3) 

— btnj — c ♦*♦♦♦♦♦;****“*♦*♦*****♦******♦********************************* 

obo C READ CHARGE DENSITY. DATA . . ■ • 

bob C QG= GAP DOPING»POSlTIVE IF SAME POLARITY AS QSS. 

000 C QP= UNDERPAD implanted doping, positive if same POLARITY AS QSS» 

""0”bcr c QSS= SURFACE state charge, CM-2 

000 C QU= UNIFOR M IMPLANTE D LAYER DENSITY >CM^2 

000' c" " CSUB=‘ SUBSTRATE DoP^^^G»CM-3 

boo READ(5»4) QSS»CSUB,QG»QP,QU . - , ^ ■ 

bob 4 FORMAT(5E10.3) 

000 r*********^:fc:»***********#******»*»*-4,**»*********** ********************* 

■ ‘bob C Esign* is used to preclude space charge INSTABILITY IN RELAXATION 

000 C IN SUBSTRAT E. IF *U » HAS WRONG SlGNr IS SET TO ZERO. 

■"'boo “ SIGN=-U0 

000 IF(CSUB*GT.0.0) SIGN = l.p ^ ■ ■ , . ^ 

000 C 3(t*),:#****#*****)(t***#***4;<:^*it:#**** + :t:*************************** 

000 C ♦3>t.»ic*********»***** **************************** ************** 

bob c CbMTuTE‘"CHARGEl)^ 

000 C ±:k*±*±**********************^1,******^********** ************** 

bob 7 “"Q=rr6E-19 

000 E0=8.85 E"*14 ■ , 

000 VSS=Q*QSS/E0*1.E“4 

000 VSUB=rQ*CS0B/EO*l.E“8 

00 b “ G0N=2.*11.7E+8*E07(Q*C^B) 
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— 00 Q_ 
opb 

OOP 

OOP 


C01=3,9E+4*EO/TOX1 

_j;i 02 = 3 - 9 E:Hf*E 0 /^ 0 X 2 

♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦**********^# + ******** 4 t**.**#: 4 c^************* 

_REMl_ELECJRQDE_yOi.TA.GES p NUMRE p nF iTERflTTONS* AMPi LOOP PRImTOUt 

* + + ** + ************************ + **stc* + **5jc***** + **;*******3*t***^** 


agjM 


vn 


irajaun: 



FORMAT(2F10.3»2I10) 

»»*»*»******»»»*» ♦**♦.* ***J*f!fe* * »»»***ifc*»»**»»» »»»»»»»»»♦»»»♦ 

DETERMINE LATERAL CONTROL INDICES 

■ ON 


■ 5 | 3 |iM«K[aHi 




(^♦♦♦♦♦♦♦*:t:****3(t*****IMPORTANTCCC IMPOPTANTCECC ♦***:*c*:»::j£**********+3t! 

C»»***»*»**:»: LG MUST BF AtvL EVFM Nl I MRER 

C LD0PE= number' OF STEPS ACROSS THE IMPLANTED REGION 

OF STEP 


C THE electrode^ POSUlVE IS TO THE RIGHT. 

Xl= NUMBER OF STpPS ACROSS THE THIN 


C LSP0X= NUMBER OF STEPS THE THiN OXIDE STARTS FROM THE LEFT EDfiE 
JL OF TH E ele ctr ode. POSITIVE TO THE RIGHT. 

REAOJ 18_) LP_i LG » LD0 PE»L SPAC> LjOKi >LSPOX 

18 FORMAT(6I10) 

I1=LG/2+1 


I2=I1+LP 

I3=I24L6 


l4=I3tLP 

L5=IMLG/2_ 

IA=I1+LSPAC 

IBFIA+LDOPE 


IC=(I2+I3)/2 

ID=I3+LSPAC 


IE=ID+LD0PE 

J 1F=I1+LS P0X 

“‘I6=IF+L0X1 ■■ 

IH=I3+LSP0X 


lK=IHtLOXl 

IL1=IF 

IF(II.'GT.IF) IL1 = I1 ■ ■ ' '■ ■ ■ 

IR1=I2 

IF<I2.GT.IG) IR1=IG 

iL2=lH 


IF(I3,6T,IH) IL2=I3 
I'R2=I4 

" lFa4;GT.IK) Tr' 2=IK 
DO 333 I=2rI5 

■ QS ( I ) =QG* ( F (T/i 1 V i r+TTi 2 » 13 » I ) +F ( 14 » 1 5 , 1 ) ) 

■ *4QP*(F.(IA»IBf I)+F.(lD>lEn) )4Qu • ' ■ 


333 CONTINUE 

QS(1)=QG . . 


/ 2 . 

/ 2 . 


QS(Il)=QS(Ii)-QG/2. 

(I2)=QS(I2) rOG/2 , 

QS(I3)=QS(I3)-QG/2. 
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“inn) — 
000 

c 


QS <14 ) 

■“OOO" 

000 

c 

c 


FIND THE EXTRINSIC DETbYE LENGTH 
**♦♦♦*♦♦*♦**♦♦***♦***♦♦**♦♦♦*♦*♦*♦***♦♦*♦***************♦♦*** 

000 

000 



OL=VTif 1 1 . TW/O/CSUS 
EDL.=SQRT(DL)*1.E4 

000 

OOO 

c 

c 


♦♦*♦***=»=♦**♦**♦***♦***♦♦*****♦♦=*♦***♦******♦***************♦* 

DETERMINE NORMAL CONTROL INDICES 

000 

000 

c 

r 


AND SET THE Y-6RID SPACING 
L1=NUMBER OF UNIFORM STEPS THROUGH SEMICONDUCTOR. 

000 

000 

c 


L2= ♦ » EXPONENTIAL * * * 

L3= ♦ » uniform steps through thin oxide 

OOO 

boo 

c 

r 


X'4= t 1 ♦ » ELECTRODE REGION. 

L5= ♦ ♦ exponential STEPS IN THE AIRSPACE REGION, 

000 

000 

c 


READ(5»302) L1»L3»L4»L5 

0.00 

000 


302“ 

FORMAT (41 10) 
Ji=Ll+l 

000 

OOO 



K1=L1 

HO=3,+TOXl/L3 . 

OOO 

OOO 



IF (HO, GT.EDU H0=E0L 

IF<HDEBY,LT.2.) HDEBY=2.0 .. . 

" boo“ 

OOO 



Y0=HDE8Y*EDL 

FO=YO/HO 

oM“ 

OOO 



Call yspac eb^eoducibility op the 

U2=J1 +L2 l^AOE. m POOR.- 

OOO 

000 



DO 10 JJ=1#L2 

OOO 

OOO 



U=J2-jJ 

H(J)= HO*ALPHA**JJ 

-OOO* 

OOO 


"IT" 

CONTINUE 

J3=U2+L3 

— ODD 
OOO 



K3= J“3~i 
DO 12 J=J2»K3 

— ow 
000 



H0X=T0Xl/O 

H(J)=HOX 

ooO 

000 


-TF~ 

CONTINUE 

I TOX= C T0X2-T0X1 ) /HOX 

000 

OOO 



I F C I fox . GT . 0 ) HT= {TeX2-f 6X1 > /I TOX 
J4=J3+ITOX 

000 

000 



K4=J4-1 

if(itox.eq.o) go to 

000 

000 



DO 21 J=J3»K4 

H(J)=HT . - . . ~ 

OOO 

OOO 


~JT~ 

CONTINUE 

J5=J4+L4 , . _ 

000 

000 


14 

K5=J5-1 
00 15 J=J4fK5 

. 000 
OOO 


15 

H(J)=TEL/L4 

CONTINUE - - - - 

000 

OOO 



"'U6=J5+L5 

K6=J6"1 - - 

000 

000 



DO 16 vJ=U5rK6 

JU~vJ“j5 , , , ■ 1 ■ ‘ ^ — j- 

000 

OOO 


16 

H ( J) = ( TEL/L4 ) * ( 2 , ♦*JJ) 

continue . - . 

000 

c 









op the 
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■ 000' 

000 
■oo'O' 

0.00 

"O06 WRITE (6» 126) Y 

000 126 FORmAT(F20.6) 

000 IF(LU.GT.IO) GO TO 127 

000 EXP0={Y-YBAR)**2/(4«*SI6MA**2) 

000 CHI.ON(LL)=EXP(-EXPO)*<H(J)+H(j-1))/2.0 

000 127 CONTINUE 

000 C * + + *'***********************)*c* + *!*:***)|c»*!*!*)»:***********^e**)t:******* + :^*** 

OOP c the following gives the WEIgHT Or THE IMPLANTED SURFACE CHARGE 
000 C WHICH IS TO BE USED AT POINTS BELOW THE INTERFACE. 

000 C *** »»*»»»***it=**»»»*********»*^ei«i»a;*********************************** 
QQQ SUM = 0 ."Q 

OOP DO 128 LL= lfl0 

000 SUM=SUM+CHION(LL) 

000 128 CONTINUE 

000 DO 129 LL=1»10 

000 CHION(LL)=CHION(LL)/SUm 

000 ' "129 ■ CONTINUE 

000 WRITE(6»131) J1 » J2r J3f J4» J5f J6 

"000 "131 FORMAT (// »10X r»J-VAL0'Es are: »» 6110"//)’ 

000 C **♦** + ♦♦♦****♦****♦♦♦*♦♦♦♦♦** + ♦*♦*♦♦♦♦♦♦* ♦*♦♦♦**♦**♦*♦***♦***4:* 

000 c compute JA»JB values WHICH Will give printout of the interface 
000 c potential distribution. 

OTTO c ******»#********>'* + *'********>^**i4i**^****i*;#!|t**'*'***4:***** + ***** 4 :* 

000 MT=(J2-JA)/JC 

“ 000 ■' JA=J2-MT*JC 

000 IF(JB.GT.J5) Jb=J5 

(JOTS c #******!t£*#******!t:**** + ****:)£**j^)je**:*:****:(:##*****2j(#******+****** 

000 c initialize potential distribution in THE semiconductor and 
OOO C ^IDE USING THE One dtmtnsionaL depletion layer estimate. 

000 c ♦2»:!«:*^:({******** + + ******3»t**i*;**4!;j. + *4:******************»5(C3»:******* 

■~?r0'0 do 40^I=2Ti“5 

000 irite=i+i 

000 IFd.EQ.I5) IRIT£=i2 

000 VI NT=Ud"l.J2)/4.4-ud. j2)/2.-»-U{tRITE»J2)/2 . 

000 XDS=AB5(C0N*VINT) 

000 XD=5QRT(XDS) 

000 Y=0.0 

000 J=J2 . ■ . ■ ■ . 

000 400 J=J-1 

000 IF(J.LE.l) GO to 402 

000 Y=Y+H(j) * ” 

000 YN=Y/XD . 

000 IF(YN.GE.i.O) GO To '402 ” ' 

000 U(I»J),=Ud»J2)»(l.o-YN)»*2 

.000 GO TO 400 

000 402 JJ=J2+1 

■ 0"00 fF(i.GE.ILl.AND.I.LEVlRl) G0‘ TO 404 

_ 000 jtpi I.GE.IL 2.AND .I,LE • I R2) 60 TO 404 ^ 

000 DO 403 J=JJ»K4 

000 , UdrJ)=U(I#J2) 

000 403 CONTINUE 

000 GO TO 406 

boo 404 bd’ 405 J=JJf K3 


I — u . 

DO 12^7_ LU=lfJ<2 

Jrj2-LL"‘ 

Y=Y-H(J) 








‘ f' M U(‘i ! 
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'UlJQ- 

000 


000 

000 

000 


U(IrJ)=U(IrJ2) 

405 CONTINUE : ^ 


406 CONTINUE 

r ♦.* *♦♦****♦ * ****t*******************************iic** 

C WRITE OUT THE STRUCTURE PARAMETERS. 

r ♦»»»:»♦♦*♦♦*»**♦*♦»*»»♦*»♦**♦**♦♦****♦*»***********♦*♦♦**♦**»♦♦♦♦»*»*=» 


WR I TE ( 6 » 345 ) TOX 1 , tOX2 t TEL » W r V 1 » V2 » LP » LG » LDOPE , LSPAC » LOXl » LSPOX » 
♦QS5fQG>QP > CSUB 


J J r r w^VtJ ^ ^ 

‘ FORMAT(/rl0Xr » ( f 0X1 fTdxaVfEL'i'V) = ♦ » 4F10 . 3» / » lOX, »■ ( Vl » V2 ) = ♦ r 2F10 . 3< 

*/>10X> « ( LP r L G »LD0PF> LS PA C »1 OXi ,\ SPQXl = » >6l6>/ ; 

*»10X» » (QSS»QG»QP»CSUB) = » »4E10.3//) 

WRITE(6f206) (U( I » J2) » 1=1 , ic) 


WRlTE(6f206) (U( I » J2) » I=ICf I5) 

♦ »*»»♦**♦ »***»» * **» ***» ***»**^***»***»»»» »**»»;<;»♦»*♦**♦♦*♦»*» 
READ RELAXATION PARAMETER VALUES 

it; *it4:**** *♦*♦**♦** *»* * ** * ************************** 

READ(5r408) 0 ‘m 1 f 0M2 » 0M3 rO'M4 » OmS’ 

F QRMAT(5F10.3) 


KT=1 

JL 0=2 __ 

JHI=J5+1 


- r ^ ■TTA jfii*i ^ J f .t. a: a: a a j 


START the relaxation PROCEDURE STARTING AT THE LOWEST ROW 
and working across to the right to the TOP ROW 


***:}<♦ ***♦♦♦♦*♦♦♦♦*♦♦*♦♦♦***♦*****♦*♦*********♦*****♦*♦*♦**** *^ 
LOOP=l 


CONTINUE 

DUMAX=0.0 


ABMAX=0,0 
DO 40 J=JL0fK2 


EA=11.7 

EB=11.7 ___ 

CALL COEFF 
0MEGA=0M1 


DO 38 I=2»I5 

ITEST=1 

IF(U(I» J+1) .EQ.0.0) GO TO 33 

LL= J2- J _ _ 

IF(LL.GT.1'0 j GO TO 3'2 

QIM P=QS(I)*CHION(Ll> „ __ _ 

QI j‘=VSUB*AVERO ( i TEStT+1 ,1-4 *qTq'tmP*'w7E0 
60 TO 35 . 


QIJ=VSUB*AVER0(ITEST> 

Call RELAXtOUrOMEGA) 


US=SIGH*U(I»J) 

IF(US. LT.O.O) U(I>J)=0.0 _ 

GO TO 3'8 

IF(U(I+1» J) .EQ.O.O) GO TO 


GO TO 30 

lF (U(I~ lfJ) .EQ.0.0) GO TO 38 

CONTINUE 

._call_mmoe_.. 

U(lf J)=U(I5»U) 

MAX=1000*ABMAX 


I F ( J . EQ . JLO . AND . MAX . GT . 0 ) JLO= J-1 

IF(MAX.EQ.O) JL0=J ^ 

IF(JL0.LT.2) JL0=2 





REPRODUCIBiLrrY OF TOE 
_ - ^CTiT/^-ixTAT, PAGE IS POOR 
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"UXTD 

000 

■W— 

c 

CONTINUE 

************************************************************* 

000 

000 

c 

c 

RELAX THE SIUCONrSILICON-nlOXlDE INTERFACE 
************************************************************* 

000 

000 


J=J2 

EA=11.7 

000 

000 


EB=3.9 

CALL COEFF 

000 


0MEGA=0M2 

000 


DO 48 1=2 f 15 

000 

000 

4a 

QIJ=VSUB*H(J-l)*W/2*+l.E-4*0*0SS*W/E0 
CALL RELAX (QlJf OMEGA) 

000 

boo 

48 

CONTINUE 
CALL MAXOF 

000 

000 

U(lrJ5=U{I5» J) 

r *♦***♦♦****♦*♦♦*♦*♦♦*+♦♦*******♦***♦********=*'***♦***♦*♦***♦****♦****** 

000 

000 

c 

c 

NEXT STEP GIVES PRINTOUT OF THE INTERFACE RtSlDUALb 
♦♦♦♦********************************************************* 

000 

000 

c 

IFdRES.GT.O) WRITe{ 6»697) DROWrIM 

♦♦♦♦♦a******************************************************* 

000 

000 

c 

c 

“relax" TTrE"'OX I dE region 

♦****4******************************+************************ 

000 

000 


JJ=J2+1 
DO 60 J=JU»K4 

000 

000 


EA=3,9 

EB=3,9 

000 

000 


CALL COEFF 
KSET=1 

oob 

000 


IF(ITOX.EQ.O) GO TO 49 
IF(J,GE.J3) KSET=0 

ot^ 

000 

49 

0MEGA=0M3 

U(i»J)=U(I5»J) 

000 

00.0 


OIJ=0.0 
DO 58 I=2rI5 

(TOTr 

000 


iF(KSEt.EQ.l) gO“ TO 5^2 
IF(I.GE.IL1.AND.I.LE»IR1) -GO TO 58- : 

000 

000 

52 

IF(I.GE.IR1.AND.I,lE.IR2) go to 58 
CALL RELAX (QIUfOMEGA) 

0015 

000 

58 

CONTINUE 
CALL MAXOF 

■'"uoir 

obo 

60 CONTINUE 

c ♦***♦**♦♦♦*♦******♦♦***♦* ****** 

000 

000 

c 

G 

RELAX THE OXlDE-AlR INTERFACE 

♦ + **************!tc************^***s)e*************************** 



J=J4 

KEH 


EA=3.9 

000 

000 


EB=1.0 
call COEFF. . 



QIU=0. 

DO 68 I=2»I5 



lF(I*GE.li.AND.I.LE*I2> GO TO 68 
IF(I.GE.I3.AND.I.LE.I4) GO TO 68 



0MEGA=0M4 

CALL RELAX (QI.J»OMEgA) 


68 

CONTINUE 
CALL MAXOF 

000 


U(lf J)="U(I5»U) 
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000 

c 


*******^*****************^***:ii*****>^**^******^*************** 

000 

c 


RE4AX THE ELECTROdE-AIRSPaCE RPGIOM 

000 

c 


♦ ♦******+***^*****************5*: + ***^*********************:*c*** 

000 



JJ=J4+l 

000 



DO 80 J=JJfK5 

obo 



EA=X* 

000 



EB=1. 

000 



CALL COEFF 

000 



0MEGA=0M5 

000. 



. QIJ=0. . 

000 



DO 78 I=2fl5 

000 



IF(I.GE.U.ANn.I.LF.I2) GO TO 78 

000 



IF(I.GE.I3.AN0.I.LE»IA) go to 78 

000 



CALL RELAX(QIJ»OMEgA) 

000 


78 

continue 

000 



CALL MAXOF 

000 



Udf J)=U(I5»J) 

000 


80 

CONTINUE 

000 

c 


*** + ************+******»*****^*********»**=tt*:(t***s,£***:*:##***** + 

000 

c 


RELAX THE AIR-SPACE REGION ArOVF THE FLFCTRODFR 

000 

c 


****«♦#**♦*♦+**♦* *%*******5jc4;*^***:(c****4:****^***=t:****!tt******** 

oob 



JJFJ5+1 

000 



DO 90 J=JJrJHI 

doo 



EA=i, 

000 



£B=1. 

000 



CALL COEFF 

000 



0MEGA=0M5 

000 



QIJ=0. 

000 



DO 88 1=2,15 

000 



CALL RELAX(QIJ, OMEGA) 

000 


88 

CONTINUE 

000 



CALL MAXOF 

000 


90 

CONTINUE 

000 



JHI=JHI+1 

000 



IF( JHI .GT.K6) JHI=k6 

000 



IFCIRES.GT.O) WRITe(6»212) DUmAy , IMAX , JMAX f JLO 

000 

c 


♦ ♦♦♦♦ + *******!(t**********!*:it***^***jt(**)tf!»;***j|i:***:*:***:4t****+4:****5t! 

000 

c 


OUTPUT CONTROL AND ITERATION CONTROL FOLLOWS 

000 

c 


this ITERATION CONTROL CAN BE REPLACED BY A CONDITIONAL 

obo 

c 


CONTROL ON EITHER dU Or OUMAX. 

000 

c 


*♦****♦*♦♦****♦♦♦♦****♦ ***5*:**:j:3(t:»:**************j4c**********:<c*** 

000 



LWRIT=LOOP/LAP 

000 



IF(LWRIT.EQ.KT) GO TO 881 

000 

207 

CONTINUE 

000 



LOOP=LOOP+l 

000 



IF(LP0P.LE.ITER) Go TO 800 

000 

c 


000 



The FOLLOWING READ MAY BE USEn FOR FINDING NEW POTENTIALS 

000 

c 


AFTER APPLYING PERTURBATIONS TO THE ELECTRODE VOLTAGES. 

000 

_£L 

♦ ♦ *♦»♦*#****** + ********** )(c»*Hc*;******:t[**:jC5tc*****5*( 

000 



READ(5»6) V1,V2»ITeR»LAP 

OOP 



IF(ITER.EQ.O) STOP 

000 



DO 883 J=J3rJ5 

000 



DO 882 I=IirI4 




iFfJ.GE.J4) GO TO 800 




U(_I.,J)=Vl*F(ILldRlrI)+V2*F(Il 2 , IR2 , 1 ) + ( 1 .-Ft II l.IRlrl) 

000 



+-F(IL2»IR2,I) )*U(Ir J) 



80 


— O'crcT” 
000 

880 

GO TO 8g2 

U(IrJ)=Vl*F(Il»I2»l)+V2*F{I3»l4,I)+(l.-F(Il»I2»I) 


000 

000 

882 

♦-F(I3»I4»I) )*U(I»J) 
CONTINUE 


000 

000 

883 

continue: 

L00P=1 


000 

000 


kt=i 

GO TO 800 


000 

000 

C DO 205 J=JA»JB»JC CONTROLS THF OUTPUT WRITE. PRINTOUT STARTS 

000 

000 

C 

881 

ON LINE JA AND GOES TO LINE Jr IN STEPS OF SIZE JC. 
CONTINUE . 


000 

000 

697 

WRITE(6»697) DROWriM 

FORMAT (2X f »RI = » rF5.3r 5x» * l = * r I3) 


000 

000 

212 

WR ITE { 6 .212 ) DUMAX » IMAx » UMAX * JLO 

FORMAT (30X» * RMAX=» »F1o,3»5X» » IMAX=» * I4*5X» »JMAX=» » 

I4»5X* • JLO=» 

000 

000 


*»I4) 

DO 205 J=JAfJB»JC 


000 

000 

206 

WRITE(6»208) J 
FORMATdOXf *J=» rl3) 


000 

000 


WRITE(6»206) (Ud»J)» I = lrIC) 
WRITE(6f206) (U{Ifj)» I=lCrl5) 


000 

000 

206 

205 

FO'R.MA'T ( i OF 10,2) 
CONTINUE 


000 

000 


IF(MAXU.EQ.O) GO TO~^r& 
DO 506 I=1»I5 


000 

000 


AUMAX=0.0 

LMAX=0 


000 

000 


DO 504 J=1»J2 
ABU=A8S{U(I»J) ) 


000 

000 


IF(ABU.GT.AUMAX) GO T<J 50? 
GO TO 504 


000 

000 

502 

AUMAX=ABU 

LmAX=J2~J 


15XrD 

000 

5F4 

CONTINUE 

DU(I)=LMAX 


000 

000 


JM=J2-LMAX 
ARRAY(I)=U(I* JM) 


00.0 

000 

506 

CONTINUE 

WRITE(6»508) 


000 

000 

508 

FORMAT C/d OX »» THE NUMBER OF "SfEPS UMAX OCCURS BELOW 
WRITE(6f206) (DU ( I ) . 1=1 , Ic ) 

SURFACE is:’/) 

000 

000 


WRITE{6»206> (DU ( I ) » I=IC » 15 ) 
WRITE(6»510) 


000 

000 

510 

FORMAT(/rlOX» »THE MAX POTENTIALS ARE:'/) 
WRITE(6f206) ( ARRAY ( I > » 1=1 r IC ) 


000 

000 

516 

WRITE(6r206) ( ARRAY ( I d I=IC f 15) 
KT=KT+1 

t 

000 

000 


GO TO 207 
END 


000 

000 


FUNCTION F(II»JJ»KK) 
F=0.0 


ogo 

000 


IF(KK.G£,II.AND.KK,LE.JJ) Fzl.O 

return 


000 

000 


END 

SUBROUTINE YSPAC 


000 


COMMON UdOOfllO) ,dU( 100) ,HdOO) »F0rL2» ALPHA » 
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"OTTQ 

OOP 

OOP 

PPQ 


PPP 

PPP 


uOP 
PPO 


PPP 



PPP 

PPP 

PPP' 

POP 

OOP 

OOP 

OOP 

poo 

OOP 


♦EA » EB~t A r.B » C » e » W rW».5TGN » DROW » qUwAX t ABMAX » 

.. „ *I>J>lMAX>J MAXf INI r 1 5 

ALPHA=i.25 

L2.=AUQgL(_0-*25*FJP±l ^q V/ALO GC ALP HA) 

2 Fl={ALPHA**L2"l.)/( ALPHA-1.) 

TRY=ARS(F1/F0) 

IFdRy.LT.l. 005. A~ND«TRY.GT. 0.995) GO TO 4 

F=F1- F0 

LL=L2-1 

F5=L2*(ALP HA**LL) /XA L PjdA.-i.iJ 

F4= ( ALPHA**L2-1 . ) / { ( ALPHA-1 . ) **2 ) 

FP=F3-F4 

ALPHA=ALPHA-F/FP 

^_T0 _2 _ 

4 CONTINUE 

■ P£TJUBN_ . 

END 

SUBROUTINE COEFF 

COMMON UdOOfilO) ,dUUOO) ,hUOO) rF0»L2 , alphas 

XE A r EBjlA » B » C r^W » VT ,_S I Gn r PR OW > h UmA X > AB M AX » 

♦i . J.IMAX. JMAXrlM.is 

A=EA*W/H(J -1) . . 

B=(EA*H( j-i)+EB*H(j) >/(2.*W) 

G=EB*W/H( J) 

C— A+2 .^B+G 

_ReJUrn__ _ 

END 

FUN CTION A VER O(ITEsT) 

COMMON U(100>110) rDU(rpo) , h ( IpO )' »F0fL2» ALPHA. 

*EA.EB.A.BrC.6rW>VT . SIGn.DRQW r pU.'^/jAX . ABMAXr 

♦I.JrlMAX.JMAX.IHf I5 
S=0.0 

~ui=ABs ( u Ti » j ) +u ( i » j-iT) ■ ■ 

IF(Ul.GT.O.O) S=0.5*H(J-1) 

iFdTESf .GT.2) GO' TO 4 ‘ - 

U2=AB5Cg(I» J)+U<I»U+1) ) 

IF(U2.GT,0,0) S=S+p.5*H(J) 

4 AVER0=S*W 

return “ 

END 

SUBROUTINE RELAX ( QlJ '.'OMEGA ) 

COMMON O(lOO.llO) .pU^OO) .HdOO) .F0.L2.ALPHA. . 

*EA r EB » a » B » C r G » W » VT r SIGN »DROW » pUMAX » ABMAX r 
*1. J.IMAX. JMAX.IM.I5 
UOLD=U(i.J) 

IRI TE=I-H 

IF<I.EQ.iS) IRITE=2 

UTIL= ( A*U ( I . JtI ) ^■B»U ( I-l . J) +B*U (IRITE. J ) -»-G*U C I .J-t-l ) +QI J) /C 
U ( I . J ) = ( 1 . -OMEGA ) *U0LD+0MEGA*UTIL 

2 DU (I )=U( I » J1-U0LD_ 

RETURN 

E ND 

SUBROUTINE MAXOF 

COMMON UUOO.nO) .nU(lpO) fH(lnO) .FO.L2. ALPHA. 

*EA. EB. A. B.C.G.W.VT. SIGN. DROW. dUmAX. ABMAX. 

*L.JZ-IiiAX_fJMAX_LlM.I5 - 

bAB=0.d 
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BEf^iaODtfCIBILITY OF THE 


CTDIT” 

000 

drow=oTO 
DO 4 I=2»I5 


IS POOR 


000 

000 

ABDU=ABS(DU(I) ) 
IF(ABDU.GT.DAB) 60 TO 2 



000 

000 

GO TO 4 
2 DAB=ABDU 




000 

DROW=DU ( I ) 




000 

IM=I 




000 

4 CONTINUE 




000 

if(dab.gt.abmax) go to 6 



000 

GO TO 8 




000 

6 ABMAX=DAB 




000 

DUMAX=DROW 




000 

IMAX=IM 




000 

JMAX=J 




000 

8 RETURN 




000 

END 




000 

qxqt 




000 

15 

80 5 

1 


000 

0.100 0.500 

0.800 

1*00 


000 

2.00E+11 3.80E-H5 

1.3SE+12 4.70E+11 


000 

1.0 20. 

240 

120 


000 

16 

8 8 

9 8 9 


000 

25 

4 10 

10 


000 

1.8 1.8 

1.5 

1.5 1.5' 


000 





000 

OEf'lD 




000 

QPCHrS TEMP.PCH 





X.DECK 

■25^-^0T7'2I^ 1 1T50 
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APPENDIX C 

FOUR-ELECTRODE PROG^ 

See Section 5, Figure 5.1 for diagram of structure. See end of 
Section 5 for description of data cards. End of listing gives 
example. 



84 


X.DL'CKr ,MINE ' “ 

B6 9 Q i/28-llt49 :03 

OOP G RUN Cj[LQ2r„EEpO ItU » JJM , » / 50 


OOP ~VaS 6»T TEMP”fF2 ’ ' 

OOP Q£LT»LIB TEMP,PCH 

OOP C **«***»****************l»:#**:»:* + **)*t3|C****5t:*j»;****3|e**:|t;*:*+)it****!(c******* 

POO C PROGRAM NO. 3 

OOP C ELECTROSTATIC ANAL'^SIS OF A CHARGE COUPLED STRUCTURE 

OOP c ♦♦***********#*****+********#*********!tr***#»3|estt***:*3»c*3|t**********!|!* 

OOP cc ■ FOUR electrode ‘cell 

OOP c itype=o : sapphire substrate 

POO C ITYPE GT, 0 J SILICON SUBSTRATE 

OOP C ******#+*****5tt***********************’*'******************^******** 

OOP c DIMENSION ARRAYS 

OOP C ^ »♦♦♦»»*★ »»♦♦ »♦»»♦» ik* »»♦»♦*»*» » »»♦»*»♦»*»♦♦♦♦»»»♦♦♦♦»»♦*♦ 

OOP COMMON U<l00rlOO)\D‘UTl00) rHa'oO) r A»B»C; s»W»EA»E8»DR0W.DUMAX»QSUB 

OOP *>ABMAXfl> J»IMAX> JMaX»IM»I5 

OOP DIMENSION DUMUPO) rOIMPUoo) 

OOP C ♦♦♦♦♦♦♦♦********4:*********)t!****:>!*****»***3t:*********^#*****!t!****** 

OOP C LGT'AND LP2“MUSf BE tVEN NUMBERS. 

POO C L61 AND LPl : GAP AND ELECTRODE WIDTH IN UNITS OF WfBOTTOM, 

OOP C LS'2 AND LP2: “ *“ »' r" » # .TOP. 

,000 C '»**♦♦** 

OOP REA0(5»2) LPl. LGi7l^"rLG2» ITYPE. MAXU 

_OOQ 2 Format (6U0) 

OOP C *♦♦♦*♦*♦*****♦♦****♦♦**♦ j(!:***>»:^*i3t:**)it**!<t***********5(ei*i^:»***i**#*#* 

OOP c LI thr ough l8 : number of steps vertically Through each of 

OOP C THE regions OeFINEd'ON DIAGRAM OF STRUCTURE. 

OOP C t******************^'***********^**^**************:)^*******#******* 

0P READ(5»3) JA.JB^JC ■ " 

OOP 3 FORMAT(3HO) 

OOP READ<5'f4> LI »r2rL3,LTfrL5'.L6»L7f L8 

__ 90 P_ FORMAT < ‘♦UP ) 

"Pop C ♦♦**=t5**»***‘***'***^**s»‘»**!»:****+*)t!******^**;(t*»*****;** + *****^*it;***#* 

POP C TOXl*TELl»TOX2»TELgfTOX3.TSI»w: THICKNESSES OF THE OXIDE> 

OOP C ELECTRODE. AND SILICON REGIONS IM MICRONS. ' ' “ ' 

_ OOP C *****^*t^* * ****^* *^********i*:^iH*^iti*ii:*^******************************* 

OOP ■ ' READ ( 5 . 6) TOXl , TEL iTToxeTteL^ » T0X3 . TSITw “ 

„ OOP 6 F O RMATC4F10.5) _ 

OOP C ***** + ******** + *********!(e+*3l>!*;. *******♦♦*♦**♦*♦♦♦♦♦♦♦*♦*!>:)♦! ♦*♦♦*♦♦♦ 

OOP c 0SS3..0S52t SURFACE STATE CHARGE AT LOWER AND UPPER INTERFACES >CMt2 
OOP C QUN: UNIFORM ImPLAnTED ION DEnSTTY»CM-2 . 

^PP- C_ CSUBt SUB STRATE DOPI NG. P OS ITIVE FOR P-TYPE.CM~3. 

OOP C ♦♦♦♦♦♦★♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦**********)|t**)((********5t:)K**3»:**i**)»:^**#* 

„000 REApJ5i8)_QSSl.QSS2»QUN.CSUB 

OOP 8 format (4E1 0.3) " ““ ' " “ 

OOP C *****»*******^i»******»**»»»i»t»;t:**=»******************************»* 

OOP C VGl»VGa»VG5»VG4 : ELECTRODE VOLTAGES VOLTS. 

. _MP. __C __LTEB?_I0TAL NO. _ OF . ITER ATIONS F OR S OLUTION. 

OOP C LAP? PRINTOUT" ON fTER." NO. LoOP/LAP = INTEGER, 

OOP C ♦»»»*»***»****» »*!» »*****»»» »**»*»4=****-»*******i>e**»***»»**»*»Ht3>*»* 

OOP REAOCSrlO) V61 .V62".VG3.VG4. ITeR.LAP " 

OOP 10 FQRMAT(4F10.3.2I16). 

000 c ♦♦♦‘♦♦♦♦*^*****:*t»*******!tt4:^:*5|:*;^***:)t)Jc*^:|c»^*:(:*+**4::,(:**:|c****:tt*^5(t***^*;*:*; 

OOP C OM N IS T HE RELAXATION PARAM ETE R IN REGION N. 

OOP C YBAR AND SiGMA ARE" THE PARAMETERS FOR THE I IMPLANTED DISTrIS. 
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000 

OOP 


000 

_0M_. 

000 

OQ-0— 


RFAD ( 5> 11) _ QM1 >O M2> 0M3_,J)Mjtj XB aR ».SI.GMA 


11 FORMAT(6F10.3) 
rn=fl.ft5t»F r>14 


Q=1.6E-19 

SieN=ltO 


000 


IF<CSUB.LT.0.) SIGn=-1.0 


000 

I2=I1+LP1 

000 

I3=i2+LG1 

000 

I4=I3+LP1 

obo 

I5=T4+LGl/2 

000 

I6=LP2/2+T 

000 

I7=I6+LG2 


000 

000 


000 

000 


I8=I7+LP2 
I9=I8-f.LG2 


IC=(l2+I3)/2 
Jl=Lt-H 


000 

000 


000 

000 


000 

000 


000 

000 


000 

do'o 


000 

000 


000 

000 


000 

000 


000 

000 


Kl=Jl-l 

JJ1=J1+1 


J2=J1+L2 
K2=J2-1 . 


JJ2=J2+1 
vJ5— J2~t‘L5 


K3=J3-1 

JJ3=J3+1 


REPRODUCIBILrry OF THE 
ORIGINAL PAGE IS -. P - QQR — 


J4=vJ3+L4 

K 4=J4~1 

jj4=J4+l 

J5=J4+L5 


K5=J5-1 

JJ5=_J5+1 

J6=J5+L6 

K6=J6"1 

Jj6=J6+l 

J7=J6+L7 


000 

000 


000 

000 


000 

000 


K7=J7-1 

JJ7^J7_+1 

VS=J7+L8 

K8=U8-1 


YSUB=I5*W 

IF<ITYPE»EQ«0) GO TO 13 


000 

odq_ 

odb 

ooj_ 

obo 

000 


ABVM=10. 
ABV1= ABS iy G1 ) 


ABV2:sABS(VG2) 

ABV3=ABS(V63) 


ABV4=ABS(VG4) 
IF(ABV1«GT.ABVM) abvw=abvi 


000 

_00Q_ 

obo 

000 


000 

000 


IF(ABV2.GT.ABVM) ABVM=ABV2 
IF(A BV3« 6T. ABVM ) . ABY^.t. AB y.3„ 


IF(ABV4.GT.ABVM) ArVM=ABV4 

VMAX=RIGN*AB VM 

C S I = 1 1 . 7E+4*e6/TS I 

CO=3 . 9E+4*E0/ ( tOXl4-TELl+TOX2 ) 


000 

000 


000 


CTOT=CO*CSI/(CO+C$l) 
QTOT= (QU N-»-QS Sl-l» E-4*C SUB*TS 1 1 *a. 


VA=VMAX+QTOT/CTOT 
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— OTTO 
OOO 

■ do‘ 6 “ “ 

OOO 

, 

Ve=Q+CSUB*ri .■7*E07TCTOT**2 ) " ■ 

PSI=V A+V 8 -SQRT { 2 . ♦ V A* Vs+VB**2 ) 

x5s=2.>j'11.7*eo*psi/(Q*csub) " 

XD=1.E+4*SQRT{A8S(X0S) ) 

OOO 

OOO 

13 

YSU8=2.*XD 

CONTINUE 



00 12 J=J2»K3 



H(J)=T0X1/L3 

■iM 

12 

CONTINUE 
DO 14 U=U3rK4 


14 

H(J)=TEL1/L4 

CONTINUE 



DO 16 J=U4»K5 
H(J)=T0X2/L5 

000 

000 

16 

CONTINUE 
DO 18 J=J5#K6 

000 

006 

18 

H(J)=TEL2/L6 

CONTINUE 

ppjl 


DO' 20 J=J6»K7 
H<vJ)=T0X3/L7 

■Ib 

20 

continue 
DO 22 J=J7»K8 

bb 

22 

H ( J ) = 0 . 2 5*H'( K7 )'*i *5**^ J~ J7 ) 
CONTINUE 

mMM 

t 

DO 24 J=Ji»K2 
H(J)=TSI/L2 

Bh 

24 

CONTINUE 

HSAP=TSI/L2*12.3/11.7 

Bm 


ALPHA=1.25 

FO=YSUB/HSAP 


26 

Fl= ( ALPHA**L1-1 . ) / ( ALPHA- 1 , ) 

TRY=ABSCF1/F0) 

000 

odo 


IF(TRY.LTa.O05.ANO«^RY.6T, .995) GO fO 28 
FA=F1-F0 

■ 'OTTO' 
000 


LL=L1-1 

F3=L1*ALPHA**LL/ (ALPHA-1 , ) 



F4= { ALPHA»*L1-1 . ) / < ( ALPHA- 1 . > **? ) 
FP=f3-F4 

000 


alpha=alpha-fa/fp 


OOO eo TO 26 

dUO 2.8“ CONTINUE 

000 DO 30 LL=1»K1 

000 J=J1-LL 

050 ___ H { J > =HS A P* ALPN A**LL _ 

000 30" CONTINUE 

000 WRl tE(6>32) _ _ 

OOO “32 FORMAT (/» 10X» » the Y-VALUES ARE:'»/) 

006 Y=0.0 ,/ ^ 

000 DO 36 J=U2rK8 

-0^0 _ WR I TE(6> 34) Y ' " 

000“ ' "“34 FORMAT<rOXVFi‘0.3) ~ 

000 Y=Y+H(J) __ 

000 DUM(jf=Y 

000 36 CONTINUE - , , 

000 YrO.O 

0^00 DO 38 LL=1^ L2 „ 

JZJ 2 -LL ■ " ' 


000 
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ow mrra:6yw)"^ 


000 
000 
000 
obb 
000 
000 
0 00 _ 
boo 
000 
boo 
boo 
000 
_o_oo_ 
000 
0 00 
000 
000 
000 
OOP 
000 
0 00 
000 
OOP 
000 

000 

000 

000 

oTo 

000 

000 

_ 0 ^_ 

boo 

000 

000 

000 

000 

000 

bod" 

OOP 

boo "■ 
000 
000 
00 ^ 
boo " 
opo 
boo' 
OOP 
000 
ooo_ 
b'b'o “ 

O OP 

boo 

OOP 

000 

OOP 

000 


Y=Y-H(J) 

DUM(J)=Y 

-QIMP(LL)=D.ISF(.Y»YbaR»Sigma)*{h(J)+K(J-U)/2». 

38 CONTINUE 


WRITE(6»34) Y 
_J=J1“UL __ 

Y=Y-H( J) 

DUM_(J)=Y 

40 CONTINUE 

SUM=0.0 

DO 42 Lt=lfL2 

SUM=:SnM+QIM P (l L) 

42 CONTINUE 

D_0_4.4J-L_=l7.L2_ _ 

Q1MP(LL)=QUN*QIMP(lL>/SUM 

44 CONTINUE 

DO 48 J=J1»J7 

D O 46 I = lrI5 _ __ 

U(I»J)=VG3*(Fdrl6»I)+F(i9/l5ri)')+VG4*F(I7»i8rI> ' 

46 Cpf£LlMUE„ _ 

48 CONTINUE 

DO 52 J=J1>K5 

DO 50 1=1 r 15 

U<I» J )=VG1*F(I1>I2>I)+VG2»F(I.^>t4»I) _ 

50 CONTINUE 
52__ CONTINUE 

fFd’fYPE.GT .0) GO TO 57 
BETA=6.28/(I5*W) 

DO 56 J=2»K1 “ ' 

X=DUM(_J) 

EXPb=EXP(BEtA*Y) ----- 

DO 1=1 d5 

u(i»j)=u(if ji)*Expo 

54 continue 
56 CONTINUE 
57_ CONTI_NUE 

WRITE (’6r5b0)’" Jlf J2»U3,"j47’j5~»’j6f J7rJ8 
500 FORMAT(/»10X» » (Jlf J2r J3» J4,J5, j7»J8) = d8I5) 

WRITE (6 » 502) LPl rLGr»LP2»L62 ' 

502 FORMAT(/>10X> » (LPl^LGl fUP2»LG2) = * r4I5/) 

WR ITE ( 6 » 504 ) TOXl r TELI r T0X2 » TeL? » TOX 3 r TSI r W 
504 FORMA T(/dOX> ♦ (TOXl r TEL2 » T0X2, TFL2» TOX3» TSI » W) = » » 7F6. 3/) 
WRITE(6r50“6) QSSl »qSS'2,qUn»CSUB 
506 FO_RMAT(/»ip_X» MOSSl »QSS2»QUNrCSUB) = * »4E10.3/) 

WRITE(6»508) Vg1»Vg2»Vg3»VG4 " 

508 FORmAT(/> 10X»< (VGl » VG2> VG3> Vf^4) = * >4FlQ ,5>/) 

WRITE (6 >510) iter# lAP»YBAR, SIGMA 

510 ..EPRMAT(/>10_X» • (ITER»LAP,YBAR»SIGMA.).=1»2I8l2F10^3jV) 

WRITE(6f512) OM1»OM2 >Om3>OM4 

. 512 .-.^.ORM AI ( OMl > QM2.>.QM:ij,.0M.4-) = Lf -4Fli}3 3 e / >_ 

MSTEP=( J2“JA)/2 

JA=J2r:MM£P*JC: ^ : . 

IF(JA.LT.I) JA=1 

. 1£JJB..*GT^J8.).. JB=U8 

C ♦♦♦♦♦♦♦***************=»*******t******************************^,+=t=* 


W TO., 

Mai PA(3fi-is--pg^^- 



1000 



JLO=kT"-3 

JHI=J 7-»1 ■ ^ . 

LOOP=i 

KT=1 — — 

CONTINUE 

JL0l=JLP:^1_ .... . 

JHI=JHl+i 

IF(J L0.LT .2) JL0=2 ; 

IFCJHI.GT.kS) JHI=K8 

DUMAX=0. 

ABMAX=0. 

IF{ITYP£.EQ»0) GO TO 660 

DO 700 J=JL0»J1 

EA=11.7 _ . . — 

”EB=11.7' 

0ME6A=0M1 

CALL COEFF 

DO 680 I=2_^I5 

"CALL aVERO 

QIJ=-U^8*Q*QSUB*cSUB/E0 

'CALL RELAX (QlUf OMEGA) 

UTEST=SIGN*U(1 > U > 

IF(UTEST.LT.O.)' U(I»J>=0^ 

co ntinue _ 

CALL MAXOF 

U(1»J)=U(I5»J) 

CONTINUE 

60 TO 720 

bo 70 J=sJLO»Kl 

EA=12._3 

■EB=12.3 

6mega=omi _ _ _ - - 

"CaLL"CO'EFF‘ 

QIJ=0t ^ 

do 68 1=2 

CALL R^AX(QIJ»0MEGAJ 

"CONTINUE " " " ‘ • 

CAU- M_A^_OF ___ 

'U(lfU')=U(l5»J) 

CONTINUE 

_RELaX_SAPPHIRE-"SILICON inter fac e _ 

_J=J1 . _ 

LL=U-J2 

EA=12,3 • , . 

' EB=U.7 

0M£ GA =0M 1 . - . I ^ 

CALL COEFF ' 

01 J =1.E-4*Q*(qIMP ( LL) /2v1.E-4*CSU B*H ( J)->-G SS2)* W/E0 

' bb’’72 I=2'fI5 '■ 

CALL relax (QIJ»0MEgA) ^ ; 

CONTINUE 

cajj- ma_xof_... : 

U(l»Jl)=U(l5f Ji) 
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.. J?EjLAX._TJiE IMPLANTED SILIC ON REGION.. 

:»:******»*********************^:t:+*5(c***:4t:)t*******!|t**’«t**************:( 

,DO.:^_jJ=5JUljdl2. . 

EA=11.7 
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X 

c 


OMEGA=:OM2 

UL=a2rJ„. 

QIJ=1.E-4*Q*(QIMP(LL)-1.E-4*CsUb*(H(J)+H(J-1))/2.)*W/E0 

_DQ.J74...I^I 5_ 

CALL RELAX (QIJ»0MEGA) 

CONTINUE 

CALL MAXOF 

. U LLt 

continue 

._ * * =L*.=t** * * *** ******* ** * ♦ ****»**★**♦=»**♦»***:< 

RELAX THE SILIC0N“Sl“0'2 INTERFACE. 

J=J2 

_ _ 

EB=3.9 

J)l^EGA=OM2 _ _ 

CALL COEFF 

LL=J-J2 

QI J=1 .E-4+Q+ ( qIMP ( lL) "1 . E-4*CsUr*H ( J-1 ) /2 . +Q5S1 ) ♦W/EO 

DO 78 I=2iI5_ . „ _ ' 

C.ALL RELAX (QI Jr OMEGA) " 

C_0MINUE_ _ 

CALL MAXOF 

Udr J2)=U(I5r J2 ) 

WRlTE(6r320) DROWriM 

FORMA TdOXr » RI = * t F8.3r5Xr t TM=* r T4 ) 

RELAX THE OX IDE R EGIO N. 

DO 90 J=JJ2rK7 . 

EA=3.9 


000 

EB=3.9 

ooo' 

Call coeff 

000 

QIJ=0. 

000 

omega=6m3 

000 

IF( J.GE. J3.i 


IF( J.GE.J5.AND.J.LE.U6) 
_ OO_80_I=2_r_I^ 

Call relax (qij .'omega') 

CONTINUE 

call MAXOF 
UC1*J)=U C I 5 r J ) 


GO TO 90 
. .P0„.S!i ..I=2.r_I_5 
IFd.GE.li.AND.I.LE»i2) 
.. IF. a.. GEL..L3...AND . I . LE .JM 
CALL RELAX (QIUrOMEGA) 
CONTINUE 


CALL MAXOF 

__0.LL? J.)_=U.LI5.»U).. 

GO to 90 


GO TO- 8? 
GO TO 86 


laEIPRODUCIBILnY OF THE 
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Go TO 84 
GO „TQ ..J8.4.. 


90 


— ^oiTQ ~SB nr=i6+i 

OOO IR=J9-1_ _ . . . 

OOQ DO 88 i=iU»IR 

000 IFd.GE.I7. A_NQ..I . UE AO. TO 8 a - 

“ 000 ' "■ Call RELAX (QIJ» omega) 

000 88 CONTINUE ^ 

000 CALL MAXOF 

0£0 ug> J)=U( l5r J ) . _ 

000 90 CONTINUE 

000 C ♦♦♦♦♦♦♦♦♦♦*******+^************«**t**Ji'*.^^^ 

~00"0 c ' RELAX 'oxide-air’ INtERFACE 

000 C ******************************* ********************************** 

000 J=J7 

__ 000 _EA=3,9 __ . 

000 e"b=i.o' 

000 QIJ=0. ■ . 

"■'ooo omega=6m3 

00 0 DO 92 I=2»I5 

000 CALL RELAX (QIJ»0MEGA) 

00_0 92_ CONTINUE _ 

■ 'ooo " CALL MAXOF 

000 __ Ud »J7)=U(I5. J7) _ _ 

boo” c ” ♦♦♦♦♦♦♦♦♦"♦♦■♦**’**********:4^*************************'*‘************** 

ooo C RELAX THE AIR-SPACe REGION. 

000 C ♦**♦♦**♦♦*♦*♦*♦*♦**♦♦*♦♦♦♦*♦:♦*♦**♦*******♦**♦♦♦***♦*♦***♦♦♦***** + 

ooo DO 96 J=JJ7»K8 

■”“0^6“ ■ EA=1. 

000 EB=1. _ _ _ 

000 ■ “ orjrb. 

000 0MEGA-0M4 ^ 

bUO DO 94 1=2.15 

000 __ C_ALL RELAX (QIJ. OMEGA) . . 

000” 94 ■continue”'' 

000 CALL MAXOF _ . 

t)W uci. jy=cj'd^f'j) 

ooo 96 continue 

000 IFdTYPE.GT.O) GO TO 99 

000 DO 97 1= 1.1 5 

■ 000 ■ ” U(I.r)=U(I»”2) 

OO Q 97 CONTINU E _ ^ 

ooo”'” 99 continue” 

000 WRITE(6.330) DUMAX. IMAX. UMAX ^ ■ - 

000 330 FORMAT(40X» *RMAX=» .F8.3.5X. »IMAX=* d4»5X. * JMAX=» .14) 

000 _ ; lwrit= loo p/lap 

boo IF(LWRIT.EQ.KT) GO TO 200 

000 98 LOOP=LOOP+l_ _ _ 

000 I'f'( LOOP. LE. ITER) Go TO ”l000 

000 READ(5.10) VG1.VG2.VG3.VG 4 . ITeR.LaP . 

000 IF(ITER.EQ.O) STOP 

000 WRlTE(6»508) VGl» V62 »Vg3.VG4 ^ 

000 W'RltE(6”.5l6‘) iTE'R.LA'PdBAR»"SlGMA 

_0£p DO 460 J=J5. J6 ' ... _ 

' OOb'"”’ ” D6’*46b‘”l=i»I5 

000 U(l.J)=VG3*(F(lfl6.I)-»-F(I9»l5.I) )-»-VG4*Fd7d8.I) 

000 460 CONTINUE 

000 480 CONTINUE 

000' ' DO ”520 j=”j3.J4 
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THH) DO 518 I=lrT5 


ooo 

ooo 

518" 

...UO.,_JJ =VS1*F ( 11 r 12 / I .) t\lQZ*F ( 1 3/ I4_f 1 ) 
CONTINUE 

ooo 

520 

^:.QNIIhLU£. ' 

ooo 


GO TO 1000 

ooo 



ooo 

ooo 

201 

WRITE(6»201) J 

format(/>iox» 


000 

000 


WRITE{6»202) 

WR ITE( 6>202) 


ooo 

202 

FORmAT(10F10.3) 

ooo 

206 

CONTINUE 

ooo 


KT=KT+1 

obo 


IF(MAXU.GT.O) GO 

ooo 


DO 212 I=1»I5 

ooo 


_UiyiAX.= Q-»_ 

bob 


DO 210 J=J1»J2 

ooo 



ooo 


IF(ABU.6T,UMAX) 

ooo 


.GD...IO 210. ... . .. 

ooo 

208 

UMAX=ABU 

ooo 


JMAX=J 

ooo 

210 

CONTINUE 


(UCIrJ) »I=1»IC) 
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GO TO 208 


000 


DUMai=U(I>J M AX) 




000 DU<I)=JMAX 

(^0 212 CONTINUE _ 

bob WRITE-{6»3d0) 

000 3^ _ FORMA T ( / tl 0 X_r • THE .VALUES OF UmA.X_ ARE :j t /)_ 

000 WRITE(6»202) (DU ( I ) » 1=1 » IC ) 

000 WRITE(6»202) (DU ( I ) » I=IC » 15) 


000 

OOP 

000 

_000 

ooo" 

o'Oo 


302 


WRITE(6»302) 

FORMAT(/rlOX» ’THE VALUES OF UmAX ARE:*»/) 

WRlTE(6f202) (DUM( I ) »'f=lVlC) 

WRITE(6*202) (pUM ( I ) > Iz iC » 15 ) 


GO TO 
END 


98 



FUNCTION F(II»JJ»KK) 

F=o_.g_ 

■ I F ( KK . GE . 1 1 VanD . KK". LE J J )' 


F=1,0 


0 00 
000 ■ ■■ 
000 
000 
000 
000 
000 
000 
000 
000 
0_00 _ 
obo ' 
000 
000 
000 
000 
0 00 - 
000 


RETURN 

e1md ■■■■ ■ ■ 

FUNCTION DISF(A»B»C) 

EXPO=-( ( A-B)/(2.*C) )**2 
DISF=EXP(EXP0) 

“return" "■ ■ 

END 

SUBROUTINE COEFF ■” 

COMMON U(lOOrlOO) fpUdOO) fH<lQ0) > A >B»C f G> W* EA »FR » PROW r DUMAX r QcUR 
♦»ABMAX»If J»IMAXf JM aX»IM»I5 
A=EA* W/H( J- D- 

B= (EA"*H ( J-1 ) +EB *H ( J )’)7 ( 2 . ‘i w") 

G=ER*W /H(J ) 

C=A+2,*B+G 

R ETUR N 

END 

..5UBR.0_UTINE AVER0_ . ; ^ 

COMMON UClOoMob) ,dU(1o' 0) »H(ioO) »A»B»C,G»W»EA»E8»DR0W»DUMAX»QSUB 
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WD- 
000 
000 
000 
000 

000 

000 QSUB=S*W 

000 RETURN 

000 ' ’ end 

000 raFORLrSI REL AX » RELAX _ 

boo COMMON uUOOrlOb) »dU<100) »H(1oO) » A»B»C,Gf WrEArEBrDROWf DUMAX»QSUB 

000 ♦>ABMAX»I»J>IMAX» JMaX»IM»I5 

000 UOLO=U(IrJ) 

OOP lRlTE=I-*-l , 

"000 IF(I.EQ.I5) IRITE=2 

0 O O UT IL= ( A*U ( I » J-1 ) +B*U < I -1 r J ) +B+U ( I R ITE f J ) +G*UJ I » J-t-D -t-QIJ) /C 

bob U(I»J)=(1 .-OMEGA)*uOLD+6ME6A*UTTL 

000 2 DU(I)=U(IyJ)-UOLD 

000 RETURN 

000 END 

000 ■ SUBROuti'NE MAXOF 

000 COMMON UdOOflOO) »oU(l00) »H{l00) »A»B»C»GrW»EA»EBrDROWrDUMAXfQSUB 

00 0 ♦ » ABMAX » I » J'f fMAXV JMaX' r Im" » 1 5* 

000 DAB=0,0 

UUTJ DROW=0.0 

000 DO 4 I=2»I5 

'000 ■ ■ ABDU=ABS(OU'(I) )■ 

000 IF(ABDU.6T,DA8) GO TO 2 

000 GO TO'4 

000 2 DAB=ABDU 

■OTJO DrOW=DU(I) 

000 IM=I 

000 4 Continue 

000 IF<DA8.GT. ABMAX) GO TO 6 

0015 60 TO 8 " ■ 

000 6 ABMAX=DAB 

WO OUMAX=bROW 

000 IMAX=IM 

000 ‘ ■ JMAX=J 

000 8 RETURN 


WO END 

000 qXQT 

M15 16 12 16 12 I 

000 5 60 5 

000 ■ 30 ’ 10 ■ ■ ' ~5 ' ' 5 

000 5 5 5 5 

000 ■ 0,4' 0.5 ' 0,4' "■ b'.a" 

OOP 0^8 do 0j_5 

000 2.00E+11 4.00E+11 l.OOE+15 

000 5.0 2p.0„ _ “lOj^O 2f^ ^ 

000 “■ 1.7 ■ 1.7 1.7 1.7" ' 100.0 

Olp 

000 (SEND 


OOP QPCHfS TEMP.PCH 


* » ABMAX f I » J » IMAX » J«AXTrM/rS 
5=0. 0 _ __ 

■ U1=A8S"( U (T f j ) +0 ( I » J-1 ) ) 
IF(Ul.GT. O.O) _S=0.5^JJH,) 
U2=ABS(U(If j)+U(IrJ+l) ) 
IF(U2.GT.0.0) S=S-t-0.5*H(J) 
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